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Abstract 


The present trend is to use high Mach number aircrafts. This necessitates the design and 
development of an aircraft gas turbine engine which could fulfill this requirement. 
Annular combustors have been found to be more suitable in high speed aircraft engines. 
There are many ways of simulating a combustor. Among these, one-dimensional 
simulation can play an important role. It gives an idea about the performance parameters 
of a combustor when there is a change in geometry and inlet conditions. One-dimensional 
analysis may be less accurate compared to 2D or 3D analysis, but the time taken and cost 
associated is considerably less. 

The present work is an attempt to simulate and analyze the fluid flow, combustion 
and heat transfer processes taking place in a high performance annular combustor to 
predict various performance parameters. Some important parameters are diffuser 
effectiveness, pressure recovery coefficient; combustor total pressure loss, mass flow 
split, mass flow distribution among various components; and the temperature distribution 
on the flame tube wall. In the diffuser part of the combustor, primarily ‘stream tube 
method’ has been used for the simulation of pre-diffliser and ‘mixing equation method’ 
for the mixing zone of annulus region. The air flow analysis has been carried out by 
solving a combination of continuity, momentum, energy, and state equations. Heat 
transfer analysis is performed by solving the equations for heat balance to obtain the liner 
wall temperature distribution, film cooling effectiveness, etc. 

The simulation results, obtained from the computer program developed, are 
compared with the experimental data and 3D analysis for the same inlet conditions and 
are found to be in good agreement. The output data are presented in tabular and graphical 
forms for better understanding of the effects of change in various parameters of 
combustor on its performance. It is expected that the present modeling of an annular 
combustor and the resulting computer program would be an useful tool for the design of 
high performance annular combustors. 



CHAPTER 1 

GENERAL BACKGROUND 


1.1 INTRODUCTION 

In modem research, computer simulation of processes and systems has become a 
powerful tool in that it saves time and is also economical when compared with 
experimental study. A proposed theory can be analyzed quickly using a computer and the 
cost of setting up an experimental apparatus can be postponed until optimization is 
achieved. However, it may be noted that simulation is only a step prior to 
experimentation and results obtained from simulation studies must be validated with 
experimental results to establish the reliability. Once validated, this tool can provide deep 
insight into the performance characteristics of the system. 

The aim of the work described in this thesis has been to develop a computer 
simulation model for the analysis of fluid flow, combustion, and heat transfer in an 
annular combustor with diffusers of making use of currently available analytical methods 
and correlations. The resulting program, used in conjunction with experimental data, will 
prove to be a useful tool for the design of high performance annular combustors. 

In a gas turbine combustor, the processes involved are extremely complex. Until 
the last decade, the design of a combustor relied heavily on previous experience and 
know-how. As a result of this, in combustor design extensive testing of the prototype is 
considered necessary. Selecting the best design from such testing is a task of considerable 
magnitude. Because of these difficulties, analysis by means of computer simulation has 
become quite popular in recent years. Further, computer simulation offers the following 
advantages; 

■ It serves as a tool for better understanding of the parameters involved and their effect 
on combustor performance. 

■ It systemizes Imowledge obtained through expensive engine testing. 

■ It reduces considerably the time-consuming tests by narrowing down the parameters 
that must be studied. 

■ It helps in optimizing the engine design for particular application, reducing cost and 
time. 

1.2 COMBUSTOR 

The function of the combustion chamber or combustor. Figure 1.1, is to accept the air 
from the compressor and to deliver it to the turbine at the required temperature, ideally 
with no loss of pressure. For the common open-cycle gas turbine, this implies the internal 
combustion of fuel and thus the problem of fuel preparation, mixing and burning. The 
fuel is commonly gaseous or liquid. Gaseous or liquid fuels are almost invariably 
hydrocarbons, the former usually being natural gas, mostly methane, with some small 
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proportion of heavier compounds such as propane and butane and some CO 2 and nitrogen, 
while the latter may range from highly refined gasoline through kerosene and 



Figure 1.1 General form of combustion chamber 

light diesel oil to heavy residual oils. It may be said that the combustion problem itself is 
seldom difficult, as all the fuels used for gas turbines are used in combustion apparatus of 
one type or another and have been used for many years. The difficulty arises in getting 
the combustion with low pressure loss in a size of combustor compatible with the high 
power to weight ratio or high output potentialities of the rotating elements. 

1.3 REQUIREMENT FOR A GAS TURBINE COMBUSTOR 
Requirements for a gas turbine combustion system include: 

i) Release of the fuel chemical energy in the smallest space (length and diameter) 
possible 

ii) Minimum pressure drop over the operating spectrum 

iii) Stable and efficient operation over a wide range of fuel-air ratios, altitudes, flight 
speeds, and/or powers levels 

iv) Reliability equal to or greater than the overall life of the engine 

v) Relight capability at altitudes for aircraft engines 

vi) Uniform temperature distribution at the inlet to the turbine stator (nozzles). 
However, many of these desirable properties are in competition. For example, 

both complete combustion and low pressure loss are contrary to small size. Hence, the 
design of a combustor is a compromise. 

1.4 TERMINOLOGY 

Certain conventional terminology is traditionally used in discussing gas turbine 
combustors. Refer to Figure 1.2 which is a cross-section of a generic diffusion flame 

combustor. 

The outer container of the combustor is called a casing. Proceeding m the flow 
direction, air exits the compressor and enters the diffuser. A portion of the air is captured 
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for the pnm^ zone by the snout. The hemispherical upstream end of the combustor 
portion IS called the dome. The fuel atomizer, surrounded by a swirler, is situated in the 
center of the dome. Air captured by the snout enters the primary zone through the swirler. 



Figure 1.2 Typical combustor cross-section. 

The body of the combustor is called the flame tube. Several rows of holes and 
slots penetrate the flame tube to admit air to the combustion process. Generally, there are 
some holes in the dome to cool its interior surface, typically called flare holes. Further 
downstream, primary holes admit air, a portion of which flows upstream into the 
primary zone, and the remainder flows downstream into the secondary (or intermediate) 
zone. At the downstream end of the flame tube, sometimes called the muff, a row of 
small muff cooling holes admit the last of the air mainly to cool the roots of the first 
stage turbine blades. At appropriate points along the flame tube, slots and louvers admit 
film-cooling air to control the temperature of the inner surface of the flame tube. The 
combustion volume is divided into three zones: primary, secondary, and tertiary 
(usually called dilution zone). 

If all the fuel and all the air passing through a gas turbine engine were mixed 
together, the resulting mixture would contain insufficient fuel to support combustion. 
More precisely, it would have a fiiel/air ratio lower than the lean flammability limit. 

Therefore, the primary function of a gas turbine combustor is to get a mixture 
which can sustain continuous combustion, and to maintain or stabilize this combustion 
over a wide range of operating conditions. This is achieved by controlled mixing of fuel 
and air. The function of primary zone is to initiate and stabilize combustion. Appropriate 
amount of fuel and air are admitted to this region to produce a near stoichiometric 
mixture, which is most favorable for both combustion and its stabilization. The 
aerodynamics of the primary zone must be arranged so as to recirculate hot combustion 
products continuously upstream to provide a continuous ignition source for fresh 
incoming mixture. Mechanical devices used to achieve this fimction are called flame 
holders. In a diffusion- flame combustor the primary zone must additionally serve as 
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region for fuel atomization and evaporation prior to combustion. The gases exiting the 
primary zone are partial combustion products at high-temperature. 

Additional air must be added to the primary zone effluent gases to more 
completely oxidize the fuel. This is accomplished in the secondary zone. If additional air 
were not added to the combustion, efficiency would be extremely low due to escape of 
incompletely oxidized species. 

The temperature of the gases leaving the secondary zone is still too high for the 
materials of the turbine section to withstand. Therefore, additional air is added in the 
tertiary or the dilution zone to reduce the gas temperature. 

1.5 TYPICAL COMBUSTOR ARRANGEMENTS 

There are different methods to arrange combustors on a gas turbine. Designs fall into 
three categories: Tubular, Annular, and Tubo-annular. 

1.5.1 Tubular (or Can type) Combustors 

Tubular combustors have approximately cylindrical flame tubes and casings (see Figure 
1 .3). Each flame tube is totally enclosed by its own casing, and the entire assembly is 
replicated in a multiple-combustor engine. 



Figure 1.3 Tubular or can type combustor 

‘Can type’ is a popular term for describing these combustors. Industrial and early 
aircraft applications predominantly employ tubular combustors. Older aircraft 
applications, particularly shaft power engines, commonly used this combustor 
configuration. Newer aircraft applications seldom incorporate tubular combustors 
because they do not provide the maximum combustion volume for a given annular space. 
Multiple igniters and/or interconnection tubes are required. However, within the 
individual combustor units, aerodynamic and combustion problems are minimized by 
tubular designs. 

1.5.2 Annular Combustors 

Annular combustors, as shown in the Figure 1.3, have a single annular flame tube and a 
concentric annular casing. This arrangement offers maximum utilization of available 
volume and is thus, widely used in modern aircraft applications. Problems related to non- 
uniform circumferential combustion (i.e. the propagation of combustion around the entire 
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circumference of the combustor from a small number of igniters) are minimal, but the 
aerodynamic performance and structural integrity are generally lower than the tubular 
designs. Achieving a uniform distribution of fuel around the annular space using a fixed 
number of fuel injectors is difficult. Maldistribution of fuel can result in non-uniform 
combustor outlet temperatures. 



Figure 1.4 Annular combustor configuration. 

Annular combustors are used mainly in aircraft type gas turbines where frontal 
area is important. This type of combustor is usually a straight-through flow type. The 
combustor outside radius is the same as the compressor casing, thus producing the 
streamline design. The annular combustor mentioned earlier requires less cooling air than 
the tubo-annular combustor, and so it is growing in importance for high temperature 
applications. 

1.5.3 Tubo-annular Combustors 

Tubo-annular combustors are hybrids of the previous two types. They have a number of 
cylindrical flame tubes contained in an annular casing (see Figure 1.5). Such combustors 
are also called “cannular” combustors. 



Figure 1.5 Tubo-annular or cannular combustor 
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1.6 MAIN COMBUSTOR COMPONENTS AND THEIR FUNCTIONS 

The turbine engine main combustor system consists of three principal elements: 

a) The inlet diffuser 

b) The dome, and snout (or cowl) 

c) The liner 

In addition, important subcomponents are- the fuel injector, the igniter, the 
combustor case, and primary air swirler. The term “combustion zone” is used to 
designate that portion of the main combustor within the dome and liner in which the 
burning of fuel takes place. These elements are illustrated below in Figure 1 .6. 



Figure 1.6 Main combustor components 

The purpose of inlet diffuser is to reduce the velocity of the air exiting the 
compressor and deliver it to the combustion zone as a stable, uniform flow while 
recovering as much of the dynamic pressure as possible. The inlet diffuser represents a 
design and performance compromise relative to required compactness, low pressure loss, 
and good flow uniformity. Early inlet diffuser designs were of the smooth curved wall or 
contoured wall type. Because of the wide variations in the characteristics of the flow field 
exiting the compressor, however, the curved wall diffuser can not always provide 
uniform, non separated flow at all operating conditions. This can become a critical 
problem in the short length diffusers required in many current systems. Consequently, a 
trend toward a dump or combination curved wall and dump diffuser designs is emerging. 
Although this design results in somewhat higher total pressure losses, it provides a known 
and constant point of flow separation at the dump plane which prevents stalled operation 
at all diffuser entrance conditions. The snout divides the incoming air into two streams- 
primary air and the other air flows (intermediate, dilution, and cooling air). The snout 
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streamlines the combustor dome and permits a larger diffuser divergence angle and 
reduced overall diffuser length. 

The combustor dome is designed to produce an area of high turbulence and floAV 
shear in the vicinity of the fuel nozzle to finally atomize the fuel spray and promote rapid 
fuel-air mixing. There are two basic types of combustor domes- bluff body and swirl 
stabilized. The bluff-body domes were used in early main combustors, but swirl- 
stabilized domes are used in the most modem combustors. 

The combustion process is contained by the liner. The liner also allows 
introduction of intermediate and dilution air flow and the liner’s cooling air flow. The 
liner must be designed to support forces resulting firom pressure drop and must have high 
thermal resistance capable of continuous and cyclic Wgh temperature operation. This 
requires use of high strength, high temperature oxidation-resistant materials and cooling 
air. 

Fuel injectors can be classified into four basic types according to the injection 
method utilized: pressure atomizing, air blast, vaporizing, and premix/prevaporising. 
Spark igniters are used to ignite the cold flowing fuel-air mixture in main burners. 

1.7 AIRFLOW DISTRIBUTION AND COOLING AIR 

The air flow distribution in, around, and through the main combustor resulting in the four 
basic air flow regions is illustrated in Figure 1.7. Effective control of this air distribution is 
vital to the attainment of complete combustion, stable operation, correct burner exit 
temperature profile, and acceptable liner temperatures for long life. 



Figure 1.7 Main combustor air flow distribution. 
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Primary air is the combustion air introduced through the dome or head plate of the 
combustor and through the first row of liner air holes. The air mixes with the incoming fuel 
producing locally nearly stoichiometric mixture necessary for optimum flame stabilization 
and operation. To complete the reaction process and control the high levels of CO, H or H2 
and unbumed fuel in the primary zone, intermediate air is introduced through a second row of 
liner holes. The reduced temperature and excess oxygen cause CO and H2 concentrations to 
decrease. In contemporary systems, the dilution air is introduced at the rear of the burner to 
reduce the high temperature of the combustion gases. The dilution air is used to carefully 
tailor exit temperature radial and circumferential profiles to ensure acceptable turbine 
durability and performance. This requires minimum temperature at the turbine root (where 
stresses are highest) and at the turbine tip (to protect seal materials). 

1.8 COMBUSTOR PRESSURE LOSS CONSIDERATIONS 

The overall combustion system pressure loss can be broken down into several 
components, as described below. 

1.8.1 Diffuser losses 

There are four components of diffuser loss- 

(i) APpj , Pre-diffiiser loss affected by inlet Mach number, inlet velocity profile and 
pre-diffuser geometry (area ratio, length and curvature). 

(ii) AP,,^^, , Small diffusion loss into cowl or snout 

(iii) AP<,p,and 

(iv) APjp , Dumping pressure losses into the outer and inner passages, affected by Pre- 

diffiiser exit Mach number, passage areas and cowl losses. 

An overall diffuser loss ( AP^ ) can be defined as the mass flow weighted average 

of AP^o^i , APgp and APjp added to the pre-diffiiser loss. 

1.8.2 Combustor losses 

There are also four components of combustor pressure loss- 

(i) APjome > Pressure drop from the cowl to the combustion chamber. 

(ii) AP„i , and 

(iii) AP;, , Pressure drop from the outer and inner liner passages, respectively, to the 
combustion zone. 

(iv) AP^t > Heat addition, or Rayleigh pressure loss. This loss is a function of Mach 
number within the combustion chamber and temperature rise. Heat addition 
pressure losses are generally very low compared to diffuser and liner pressure 
drops. 

An overall pressure drop ( APij^^er) ^an be defined as the mass flow weighted average 
of AP , AP, , and AP„, added to the heat addition loss. 

dome ’ U ’ 01 
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It is desirable to minimize dififiiser loss since this is a parasitic loss. On the other 
hand, liner loss is used to promote uniform air distribution and provide high swirler and 
dilution air velocities for uniform and intense mixing within the combustor. 

1.9 LINER COOLING TECHNIQUES 

Cooling air must be used to protect the burner liner and dome from the high radiative and 
convective heat loads produced within the burner. This air is normally introduced through the 
liner so that a protective blanket or film of air is formed between the combustion gases and the 
liner wall. 

The effectiveness of cooling technique is quantified by the cooling effectiveness (s), 
defined by 


T -T 

8 c 

where Tg, Tm, and Tc are the mainstream gas, average metal, and cooling air temperatures, 
respectively. 

As future combustor exit temperature requirement increases, the percent of combustor 
air available for cooling decreases requiring increased cooling effectiveness. 

1.10 PREVIOUS EFFORTS 

A number of attempts have previously been made to calculate the air flow and heat transfer in 
gas turbine combustors. Grraves [1] and Grobman [2] analyzed the pressure loss and air-flow 
distribution in tubular combustors with uniformly tapered (or constant cross-section) flame 
tubes. This work included the effects of momentum transfer between the gas streams in 
annulus and flame tube, annulus wall friction, heat release, hole discharge coefficients, and 
compressibility. Flow in diffuser was not considered, and instantaneous mixing between cold- 
air jets and the gas in the flame tube was assumed. No heat-transfer analysis was undertaken. 
The results of these calculations exhibited many of the effects observable in practical 
combustors. 

Samuel [3] used the same general approach as Grobman, with instantaneous mixing 
between cold-air jets and the main gas stream. A simple diffuser analysis was included, and 
both tubular and annular geometries could be treated. In the case of annular combustors, a 
calculation followed the three parallel streams in the inner and outer annuli and through the 
flame tube, and iterated on the initial mass-flow split between these streams until boundary 
conditions at the end of the combustor were satisfied. A similar approach is used in the 
present work. The output of Samuel’s program was used as input for the heat transfer program. 

The most comprehensive heat-transfer analysis of aircraft-type combustors that has 
been published to date is that of Lefebvre and Herbert [4]. This involved the solution of a 
heat-balance equation for each element of the wall, taking into account radiation from flame 
to wall and wall to casing, and forced convection on the inside and outside of the wall. Net 
heat transfer in the longitudinal direction was assumed to be negligible compared with radial 
heat transfer. Reasonable agreement with limited experiment was obtained. 

Another heat-transfer analysis was carried out by Tipler [5] whose work dealt with an 
industrial combustor and included the effects of radiation transfer in a longitudinal direction. 
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Another heat-transfer analysis was carried out by Tipler [5] whose work dealt with an 
industrial combustor and included the effects of radiation transfer in a longitudinal direction. 

Sovran and Klomp [6] provided a way of generalizing straight-walled two- 
dimensional and annular diffuser performance. An almost identical approach was used by 
Reneau et al. [7], who extended the work to cover different inlet boundary- layer thicknesses; 
this aspect was particularly important for gas-turbine combustor diffusers. For the first time, 
therefore, there existed a set of maps which enabled the performance of straight-walled 
diffusers to be predicted with fair certainty for any length and area ratio and inlet boundary- 
layer thickness, at least within the range of interest. 

In the field of heat transfer, Spalding [8] correlated the results of 9 experiments to 
produce a new correlation for film-cooling heat transfer. Schirmer and Quigg [9] produced the 
most useful work yet published on the effect of pressure on radiation from luminous flames. 

The works mentioned above, and others, provide useful new tools for combustor 
analysis, and open the way for removing some limitations that have previously limited the 
utility and validity of such analyses. 

1.1 1 SCOPE OF THE PRESENT WORK 

Present thesis work is devoted to carry out an overall analysis of an aircraft annular 
combustor operating under different conditions. For this purpose, an attempt has been 
made to develop a computer program for one-dimensional analysis. The computer 
program is intended to achieve the following objectives- 

To analyze the air mass flow and pressure distribution in a combustor whose 
geometry is given, for particular inlet conditions. 

> To compute the temperature distribution on the flame tube wall. 

> To find out pressure loss in various components of an annular combustor, viz., 
diffuser, annuli, flame tube, and swirler etc. 

> To analyze the mass flow split for inner annulus, outer annulus and flame tube. 

> To analyze the mass flow distribution in cooling rings and cooling holes. 

1.12 OUTLINE OF THE THESIS 

The main work of this thesis has been described in various sections under Chapter 2 through 
Chapter 6. Chapter 2 is devoted to describe the methodology for the one-dimensional analysis 
of the combustor. It discusses computer model development in which the input data 
requirement and methods of the geometrical indexing are given for calculation at various axial 
locations of the combustor. Chapter 3 deals with the analysis of combustor diffuser. In this 
chapter different methods of diffuser calculation are discussed. An outline of the calculation 
procedure is explained in a detail. Chapter 4 discusses the methods for air flow analysis in the 
annuli and flame tube. Prediction of air mass flow split in the annuli and air mass flow 
distribution through various cooling holes and cooling rings are included. Chapter 5 is 
devoted to combustor heat transfer analysis. The primary objective of heat transfer analysis is 
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to establish the axial distribution of temperature along the flame tube wall. To serve this 
purpose, the heat balance equation is solved to determine the unknown temperature. Chapter 6 
presents the overall result obtained from the computer simulation. The results are compared 
with the available experimental data and 3D analysis for validation of the model and the 
resulting program. 



CHAPTER 2 


MODELING METHODOLOGY AND PROGRAM DEVELOPMENT 


2.1 INTRODUCTION 

This chapter is an introduction to the combustor simulation and the assumptions made and the 
methodology used. It is intended to provide a frame work for the detailed discussion that 
follows in the succeeding sections. 

2.2 MAJOR ASSUMPTIONS 

The following assumptions are made for the present: 

> Flow and combustion conditions are steady: no quantity varies with time. This 
assumption is reasonable, but it precludes the treatment of transient phenomena such 
as combustion oscillation. 

> Combustors are of annular cross section. This precludes the direct simulation of tubo- 
annular (cannular) and tubular combustors. While flow conditions for these types may 
be approximated by appropriate selection of dimensions, other parts of the simulation 
(for example, radiation heat transfer) will be of reduced accuracy. 

> Conditions are uniform around the circumference of the combustor; i.e. there is no 
circumferential variation in any quantity. This assumption is a major simplification; in 
practice, cyclic variations around the circumference occur, corresponding to the fuel- 
injection points. The result of this analysis would be useful, however, in indicating 
output parameters averaged around the circumference. 

> The flow can be treated as quasi-one-dimensional. To represent the complex flow 
pattern existing in a real combustor is beyond the scope of the present worL Instead, 
the problem is broken down into two parts: flow in the main gas streams (one- 
dimensional) and flow and mixing of jets (predicted from correlated experimental 
data). 

> Radial-flow or reversed-flow in the combustors is not analyzed. 

> Flow conditions in the combustor are within the range of validity of the correlations 
used. The correlations are based on experimental data presently available. Outside 
their range of established validity, accuracy of prediction is likely to be reduced. 

2.3 INPUTS TO THE COMPUTER PROGRAM 

Input quantities required by the computer program are the following types; 

(i) Combustor geometry data 

(ii) Inlet flow conditions. 

The following information must be supplied to specify the combustor geometry: 

^ Dimensions of diffuser walls, snout, dome, flame tube, and casing 

> Hole and cooling-ring (equivalent cooling-slot ) data: 
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a) Type of holes 

b) Number of holes in a row 

c) Axial locations; inner or outer wall 

d) Area (or cooling-slot height) 

> Identification of secondary-hole row (marks end of primary zone) and first hole row or 
slot in the flame-tube wall, as distinct fi-om the dome (marks end of diffuser and start 
of combustor annulus) 

> Specification of swirler (optional): 

a) No swirler 

b) Swirler designed within program 

c) Swirler dimensions supplied as input 

2.4 GEOMETRICAL INDEXING SYSTEM 

Four interrelated geometrical indexing system are used within the program- 

1. Input stations (i.e. geometric input stations) 

2. Hole rows stations and cooling rings (equivalent slots) 

3. Diffuser stations 

4. Flame tube calculation stations 

These are illustrated by the sample numbering system on the Figure 2.1: 



Figure 2.1 Geometrical indexing system. 
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2.4.1 Input Stations 

The shapes of the walls of the diffuser, snout, dome, flame tube, and casing are described by 
their axial and radial dimensions at geometric input stations. Axial dimensions are measured 
from the diffuser inlet (compressor exit); radial dimensions from the engine center line. 
Station numbers begin with station 1 at the compressor exit. 

2.4.2 Hole Rows 

Holes are indexed in groups referred to as hole rows. A hole row comprises holes of identical 
size and shape, centered at the same axial location, and spaced at equal intervals around a 
single wall of the flame tube or dome. Up to six hole rows may be specified at a given axial 
location: hole rows may be interspersed on the same wall or paired on opposite walls of the 
flame tube. Indexing starts with the first hole row down stream of the compressor exit. 

The indexing system for slots is similar to that for holes, except that slots are 
circumferentially continuous (as mentioned above) and hence only one slot may be placed on 
a given wall at any axial location. Non-continuous slots may be approximated by continuous 
slots of equal area. Slots and holes are indexed together. 

Dummy hole rows (holes of zero diameter) may be specified as input, for purposes of 
arranging the location of calculation stations. 

2.4.3 Diffuser Stations 

The calculations in the diffuser make use of the input-station indexing system. In addition, 
labeled as follows in the above figure 2. 1 : 

1 Entrance to the diffuser (compressor exit), 
r -^Outlet to the pre-diffuser or dump diffuser inlet. 

2- > Last geometric input station upstream of the start of the snout (or cowl). 

2' -^First geometric input station on the snout. 

3- > Last geometric station where diffusion takes place. 

4- > Entrance to the annulus (diffuser exit). This station coincides with the upstream 

edge of the first hole row in the flame tube, as distinct from the dome. 

2.4.4 Flame-Tube Calculation Stations 

The primary zone of the flame tube is bounded by the dome and swirler, the flame-tube walls, 
and (at its down stream end) by a plane passing through the upstream edges of the secondary 
holes. This zone is treated as a homogeneous stirred reactor; no calculation stations are used. 

Calculations in the annulus and in the flame tube dovra stream of the primary zone are 
related to a system of calculation stations located as follows. 

1 . At the upstream edge of each hole row that is centered at a distinct axial location. Thus, 
for hole rows paired on opposite flame-tube walls only one calculation station is 
assigned. 

2. At each dummy hole row. 

3. At the axial location of each cooling slot i.e. cooling ring having a distinct axial 
location. 
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4. At up to five axial locations spaced at a specified intervals down streanm of each 
cooling-slot location. 

5 . One station at the very end of the combustor. 

Indexing of calculation stations starts with the first hole row in the flame tube (as 
distinct from the dome) down stream of the compressor. The first calculation station thus 
coincides with the end of the diffuser annuli (Diffuser Station 4) and the start of the 
combustor annuli. 

2.5 OVERALL STRUCTURE OF COMPUTER MODELING PROGRAM 

The present computer model comprises many subprograms; these can be grouped by 
function into three major sections as follows: 

1. Diffuser subprogram 

2. Air flow subprogram 

3. Heat-transfer subprogram 

2.5.1 Diffuser Subprogram: 

The diffuser subprogram receives inputs of geometry, inlet flow conditions and certain 
quantities related to the calculation procedure; it also asks a first estimate of mass flow 
among the cowl (i.e. snout) and the two annuli. 

Using these inputs, the diffuser subprogram performs the following functions; 

(a) Determines Diffuser Performance Parameters: 

> Ideal pressure recovery coefficient 

> Actual-pressure recovery coefficient 

> Effectiveness 

(b) Evaluate the mismatch at the snout or cowl 

(c) Calculates flow conditions on the dome and at the end of diffuser annuli 

Within the diffusing regions the three calculation methods may be used which are; 

1 . Stream tube method 

2. Empirical data method 

3. Mixing equation method 

These methods will be described in detail in the next chapter. 

2.5.2 Air Flow Subprogram: 

The air flow subprogram receives the following inputs: 

• Geometry of walls, holes and swirler at appropriate geometric stations. 

• Fuel data 

• Initial estimates of mass flow-split among snout and annuli 

• Static pressure and static temperature on the dome and at the end of diffuser 
annuli. 
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With the use of these inputs, the air flow subprogram performs the following function: 
Calculates the flow conditions in the primary zone and at calculation station in 
the flame tube and annuli. These flow conditions are static pressure and 
temperature, mass flow rate and velocity 

> Calculates the combustor total pressure loss 

> Directs the iterative process by which the diffuser and air flow subprogram 
together arrive at the correct flow split among the snout and two annuli. 

2.5.3 Heat-Transfer Subprogram: 

The primary objective of heat transfer subprogram is to establish the axial distribution of 
temperature along the flame tube walls. 

Geometry of the flame tube and casing, and the axial distribution of velocity and 
temperature of the flame tube gas are taken as input in the heat transfer subprogram from 
the air flow subprogram. The additional data may be specified as input to the program. 

In operation, the subprogram evaluates various heat flux components at a point on 
the flame tube wall in terms of the wall temperature; the heat flux components which are 
considered include the following: 

i) Convection from the flame tube wall gas 

ii) Convection to the annulus air 

iii) Radiation from the flame 

iv) Radiation to the outer casing. 

The heat balance equation is then solved to determine the unknown temperature. The 
operations are performed at calculation stations along the flame tube wall. 



CHAPTER 3 


COMBUSTOR DIFFUSER ANALYSIS 


3.1 FUNCTION OF THE COMBUSTOR DIFFUSER 

In the combustor inlet of an aircraft gas turbine engine, high velocity air fi'om the 
compressor flows into the diffuser, where a considerable proportion of the inlet velocity 
head is converted to static pressure before the airflow enters the main combustor. Modem 
high through-flow turbine engine compressors are highly loaded and usually have high 
exit Mach numbers. With high compressor exit Mach numbers, the velocity head at 
compressor exit may be as high as 10% of the total pressure. The function of the diffuser 
is to recover a large proportion of this energy. Otherwise, the resulting higher total 
pressure loss would result in a significant loss of power and fuel efficiency. The diffuser 
performance must also be insensitive to inlet velocity profiles and geometrical variations 
of the combustor relative to the location of the pre-diffuser exit flow path. 

Low diffuser pressure losses with high inlet Mach numbers are more readily 
achievable with increasing length. But diffuser length must be short to minimize engine 
length and weight. A good diffuser should have a well considered balance between the 
conflicting requirements for low pressure losses and short engine length. 

3.2 DIFFUSER REQUIREMENTS 

Important design requirements for combustor inlet diffusers that must accept high 
compressor exit Mach numbers are as follows; 

(i) Low pressure losses. In general diffuser pressure loss should be less than 40% 
of the compressor exit velocity head. 

(ii) Short length. Configurations that have special features, such as annular splitter 
vanes can be used, in some cases, to reduce length. 

(iii) No flow separation, except in dump regions. 

(iv) Uniform flow both circumferentially and radially. 

(v) Dynamic flow stability at all operating conditions. 

(vi) Insensitivity to changes in compressor exit flow patterns or exit flow 
conditions. 

3.3 DIFFUSER TYPES 

Diffusers are of different types based on requirements and engine operating conditions. 

The main types are: 

3.3.1 Step Diffuser 

This configuration is mechanically simple, aerodynamic^ly efficient, and has a relatively high 
tolerance to distorted inlet velocity contours and dimensional tolerances. This type 
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of diffuser is also known as annular step diffuser. This thesis deals with only this type of 
diffuser keeping in view the indigenous Kaveri gas turbine aircraft engine (proposed to 
be used in LCA) development efforts, and will be discussed in detail 

3.3.2 Controlled Diffuser 

The controlled diffuser is designed to diffuse the compressor discharge flow smoothly 
with no flow separation at any operating conditions. This design characteristic imposes 
limitations on the overall dif^ser area ratio and flow turning capability of the diffuser 
walls. This type of diffuser has a short pre-diffuser and long combustor cowling that 
extends upstream to the pre-diffuser exit plane. The controlled diffuser has low pressure 
losses in the inner and outer passages, but the center flow pressure losses are usually 
higher that those in the annular step diffuser. Advanced combustion system trends limit 
the use of controlled diffuser. 

3.3.3 Multiple Passage Diffuser 

An annular splitter vane is used in this concept to reduce the length required to obtain the 
desired diffuser area ratio and to provide high-energy flow to the combustor dome region. 
With a high area ratio pre-diffuser and low dumping losses, the overall pressure losses for 
this configuration are very low. 

3.4 ANNULAR STEP DIFFUSER 

This type of diffuser mainly has two flow regions to fulfill the requirement of short 
diffuser length and low diffuser losses. The configuration is illustrated in Figure 3.1. 



Figure 3.1 Annular step configuration. 
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The main regions are as follows; 

3.4.1 Pre-diffuser 

A large fraction of the compressor exit velocity head is recovered in the pre-diffuser 
sections of the combustor inlet diffuser. This reduction in velocity head results in reduced 
dumping losses at the end of pre-diffuser and reduced parasitic drag losses from frame 
support stmts and fuel nozzle stems. The pre-diffuser must be carefully designed to avoid 
the possibility of flow separation ahead of the step regions. If flow separation occurs 
within the pre-diffuser section, the separated regions are usually unstable and occur at 
isolated locations around the annulus near the exit end of the pre-diffiiser, sometimes at 
particular engine operating conditions. These randomly located flow separations cause 
severe circumferential distortions of the flow entering the combustor, which may result in 
hot streaks at the turbine inlet plane. These flow separations can also lead to dynamic 
instability. 

Pressure losses in the pre-diffuser are relatively low and nearly all of the losses 
are caused by momentum deficiencies in the wall boundary layers. The velocity profile 
becomes more center-peaked and the boundary layers become weaker as the flow moves 
down stream. For the same average velocity, a center peaked profile has a higher mass- 
weighted velocity head than a flat profile. With more of the energy in the velocity head, 
the static pressure recovery is reduced below the ideal levels. High turbulence levels 
produce more uniform profiles and reduce the pressure losses. 

3.4.2 Dump-diffuser (or Diffuser Step Regions) 

Annular step diffusers are designed to have backward facing steps at the end of the pre- 
diffiiser. The typical regions are illustrated in Figure 3.2. The flow at the pre-diffiiser exit 
plane separates from the sharp comers at the trailing edge of the pre-diffriser and curves 
around the combustor cowling to the combustor liner passages. At the base of the dump 
region the flow is recirculating at low velocities. This low-momentum recirculating flow 
mixes with the higher energy flow leaving the pre-diffiiser, resulting in momentum, 
mixing losses as the flow goes into the liner passages. Front of the reattachment point of 
the stagnation streamline, flow is reversed and goes into the recirculating region, and aft 
of this point, the flow goes downstream. If stable flow patterns exist in the recirculating 
region, the forward flow at the stagnation point is balanced by the flow mixes into the 
mainstream 

A considerable proportion of the pressure loss in the annular diffuser is due to the 
mixing losses in the dump region. These losses are difficult to estimate because the 
geometry of this region is complicated. The theoretical dumping losses for step increase 
in passage area can be obtained by applying the momentum equation with the assumption 
of uniform, steady incompressible flow with negligible friction. 
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where A3.1 and A 3, 2 are the passage areas at stations 3.1 and 3.2, which are located 
upstream and downstream of the step respectively as shown in the Figure 3.2 . 



Figure 3.2 Recirculating flow patterns in dump diffuser region 

One major difficulty in the analysis of dumping loss is to select an appropriate 
value for A 3 . 2 , which is the effective cross-sectional area of the passage at the unknown 
point where the flow reattaches to the wall, i.e., the stagnation point. Dumping losses can 
also be affected by inlet velocity profiles and inlet turbulence levels, and by physical 
spacing between the exit plane of the Pre-diffuser and the cowling leading edges. If this 
spacing is too small, the flow into the dump region will be accelerated around the 
cowling leading edges and the pressure losses would be higher. 

3.4.3 Dome Regions: 

Ideally, the region within the combustor cowls, upstream of the dome, should be a 
plenum chamber with very low velocities and a uniform static pressure to provide 
uniform flow into the swirl cups and other openings. The upstream dome region of the 
annular step diffuser configuration, as shown in Figure 3.3, approaches this ideal concept. 
The area ratio of the center flow stream, fi-om the exit plane of the pre-diffuser to the 
combustor dome region is very high, typically about ten to one. A large amount of 
diffusion occurs at upstream of the cowling opening as free stream diffusion. The 
remaining diffusion takes place inside of the cowling where the dome blockage effect 
spreads the flow out to the inner wall surfaces of the cowlings. 

The cowling opening must be properly sized to achieve sufficient fi'ee stream 
diffusion without flow spillage from the cowling. If the opening is too small, the static 
pressure recovery in the dome region will be reduced. But if the opening is too large, 
some of the flow entering the cowling will turn around and flow out around the upper or 
lower cowling leading edge, as shown in Figure 3.3. This flow spillage firom the cowling 
will reduce the pressure recovery in the dome region and will cause higher passages 
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losses, because the spillage flow, with low momentum, mixes with the higher velocity 
passage flow as this flow goes into the combustor liner passages. In some cases, this 
spillage flow may cause flow separation from the outer cowling surface. This depends on 
the cowling lip contours and on the contours around the outer sur&ces of the cowling. 
The cowling contours should have smoothly accelerating flow on the outer surfaces from 



Figure 3.3 Diffuser with flow spillage around the outer cowling leading edge 

the stagnation point at the leading edges of the combustor liner passages. A semi- 
elliptical shape, with the long axis horizontal approaches the ideal configuration. 

Pressure losses in the flow that goes into the dome region are very difficult to 
determine. These losses are usually low and the shape of the inlet velocity profile heavily 
influences them. 

3.4.4 Combustor Liner Passages: 

Static pressure in the combustor liner passages should be as high as possible to provide 
the uniform flow in the combustion zone. The static pressures depend upon the passage 
area ratio and the velocity profile in the passages. If the passage area is too low, the 
velocity will be high and the static pressure would be lower. However, higher passage 
area ratio cause the dumping losses to be high and the flow leaving the dump region may 
not have reattached to the casing wall before entering the liner passage. In this case, the 
velocity profile in the passage would be much distorted, which would result in low static 
pressures. To avoid dynamic instability the flow should be reattached. If the flow in one 
passage is not attached, a disturbance which momentarily increases the flow in that 
passage would tend to fill the passage, resulting in reattachment. Under this condition, the 
passage static pressure would increase with increasing flow to the passage. This is an 
unstable situation that must be avoided. 

Computer codes to analyze compressible viscous flow are continuously being 
improved and are being applied in diffuser analysis. This thesis is a bid in the same 
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direction and a computer program for one-dimensional analysis of diffuser and main 
combustor has been developed. 

3.5 ANALYTICAL METHODS AND PROGRAM DEVELOPMENT 

FOR THE DIFFUSER 

As mentioned in the Chapter 2 , the three methods for diffuser analysis are; 

3.5.1 Stream Tube Method 

In this method the flow passage is divided into N stream tubes within each of which flow 
is assumed to be uniform and isentropic. In setting up the stream tubes, the distribution of 
boundary layer displacement thickness supplied as a program input is used to adjust the 
flow area. The development of the flow (static pressure and velocity) is computed for 
each stream tube at each station along the diffusing passage. 

Boundary layer calculations are performed to obtain a revised estimate of 
boundary layer displacement thickness, which is then used in a revised stream tube 
analysis. This process is continued to convergence of the displacement thickness. The 
result of this method is the velocity distribution, static pressure and boundary layer 
blockage at the outlet from the diffusing passage. 

3.5.2 Empirical Data Method 

This method is based on the direct evaluation of the effectiveness of a particular diffusing 
passage. There are two ways this can be done- 

a) Diffuser effectiveness may be supplied directly as a program input 

b) Diffuser effectiveness may be calculated by the program fi-om correlation of 
experimentally measured effectiveness versus diffuser geometry that is contained 
in the library data. 

The re mainin g diffuser performance parameters like ideal and actual pressure 
recovery coefficients and the outlet conditions are calculated directly from the 
effectiveness, the geometry and the assumptions that the total pressure is constant. 

3.5.3 Mixing Equation Method 

The mixing- equation method is always used in the mixing regions. It takes into account 

the following effects occurring in these passages; 

i) Pressure loss due to sudden expansion or contraction at the snout or cowl. 

ii) Mixing in the diffusing passage from dome to the annulus diffusion region 

iii) Pressure loss due to curvature of the flow passage from dome to the annulus 

diffusion region . j.-rr- p 

Above three methods of diffuser analysis are separately useful m different part of 
the diffuser. Their usefulness and application in the model development is described in the 
succeeding sections of the diffuser analysis 
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3.6 OBJECTIVES OF THE DIFFUSER ANALYSIS SUBPROGRAM 

The objectives of the subprogram are: 

♦ To provide the necessary inlet flow conditions to the annuli; static pressure and 
temperature the specific properties calculated, although all others can be derived 
from these two and the mass flow. 

♦ To determine the diffuser performance parameters of effectiveness, ideal pressure- 
recovery coefficient, and actual pressure-recovery coefficient. 

♦ To determine whether separation occurs in the diffuser and, if so, at what points. 

♦ To evaluate the magnitude of the flow mismatch (for diffusers with a snout). 

3.7 MAJOR ASSUMPTIONS IN THE DIFFUSE SUBPROGRAM 

The assumptions employed in the diffuser analysis common to all calculation options, 
referring to Figure 3.4, are that; 

• The flow in the upstream part of the diffuser (l-T, r-2) is unaffected by the flow 
downstream of Station 2. 

• If there is a snout (or cowl), it is reasonably well matched to the flow; i.e., the 
streamlines in the portion 1-1', r-2 which bound the flow required by the annular 
passages should have radial locations at Station 2 nearly equal to those of the snout as 
shown in the Fig. given below- 



Figure 3.4 Flow match and mismatch at snout station 

If the flow is not well matched, it is certain that the flow in the upstream part of the 
diffuser will be affected by flow in the down stream part, in which case the results of 
the diffuser subprogram are highly questionable. 

• If there is a snout i.e. cowl, its performance as a diffusing passage is specified by an 
input value of the total-pressure-loss coefficient. This is based on the presumption that 
the flow in the snout (or cowl) will be dominated by either the colander (or orifice plate) 
or the downstream resistances supplied by the swirler and dome holes; in either case, it 
is not expected that the total-pressure loss due only to diffusion in the snout (or cowl) 
will be of sufficient magnitude to warrant a separate analysis. 

• If there is no snout (or cowl), the performance of the passages 2-4A and 2-4B is 
specified by a single input value of total-pressure-loss coefficient between Stations 2 
and 4 (i.e., ratio of total-pressure loss between stations 2 and 4 to the dynamic l^ad at 
station 2). This is based on the presumption that in such cases the resemblance of these 
passages to a diffuser will be remote. 
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• The static pressure across any diffusing passage is constant across any section normal 
to the passage centerline. This assumes that the flow will be subjected to mild 
curvatures, and it seems to yield reasonable results for well-designed diflEusers. 

3.8 DEFINITION OF TERMS USED IN T HE SUBPROGRAM 

The terms fall in two categories: flow parameters, and diffuser-performance parameters. 

3.8.1 Flow Parameters 

The flow parameters which are most frequently used are the area-average velocity, the 
mass-average velocity, the ratio of equivalent free-flow area to total flow area, or blockage, 
and the dynamic head. All of these parameters refer to the flow at a given axial location, 
where the static pressure is assumed constant across any section normal to the passage 
centerline. 


> Area-average velocity, u , , is defined by 

Ua =- judA : (3.1) 

where, u is the velocity component in the direction of the passage center-line (normal 
components are neglected) and ‘A’ is the passage cross-sectional area normal to the 
passage center line. For incompressible flow, u ^ is the velocity the flow would have if the 
profile were uniform across the entire passage. 


> Mass-average velocity, u jjj, is defined by 

f uMA 

u,= i 

m A 

1 udA 

0 


(3.2) 


For incompressible flow, it is the velocity which when multiplied by the mass flow would 
yield the total momentum flux. 


> Blockage is 1-E, where E is defined by 



0 


(3.3) 


For incompressible flow, E is the ratio of the area which would be occupied by a uniform 
flow, with a total pressure equal to the mass-averaged total pressure of the actual flow, to 
the actual area of the passage. A related term which is frequently used is the profile 



25 


parameter, P , which is merely the reciprocal of E. a derivative of the general inlet blockage, 
1-E, is the inlet boundary layer blockage, l-E^^ , where E is defined by 


1 


(3.4) 


where A^j is the area occupied by the boundary-layer displacement thickness. If the flow 
outside the boundary layer is uniform (i.e., u is constant across the passage) then E = E 


^ Dynamic head based on mass-average velocity, q , is defined by- 

q = Un> ^ ^ 

z, 

where p is the average density across the passage, mathematically 
m 


(3.5) 


P = = 


(3.6) 


UmEA 

For incompressible flow, q the actual mass-average dynamic pressure; for compressible 
flow, it is treated as a reference quantity. 

3.8.2 Performance Parameters 

The diffuser performance parameters used herein are the ideal incompressible-flow 
pressure-recovery coefficient, the actual pressure-recovery coefficient, and the diffuser 
effectiveness. 

> Ideal pressure-recovery coefficient is defined by 


Cp, 


ideal “ 


. (P2-_Pl)n 

qi 


(3.7) 


where the subscripts 1 and 2 refer to inlet and exit, respectively. 

> Ideal incompressible-flow pressure-recovery coefficient, Cpi , defined by 


r * 


Cp.i=l- 




(3.8) 


where, — is the passage area ratio. This expression in fact represents the maximum static- 
'll 

pressure rise that can be achieved in incompressible flow if the velocity profile at the 
diffuser inlet is uniform. If, however, the inlet velocity profile is non-uniform an additional 
increase in static pressure may be expected due to mixing process. 

> Ideal pressure-recovery coefficient in presence of mixing. Cm , is given by 
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C = — — 


2(p,-l)+l-fA 

I Ao 


where Pi 3 1 for uniform profile and 3i>l fijr non-uniform profile. 

The Equations (3.8) and (3.9) are more appropriate only for incompressible flow. 

> Actual pressure-recovery coefficient, Cp , is defined by 
. _P2-Pi 

Cp = 


(3.10) 


> Diffuser effectiveness (4^) ; 


The diffuser effectiveness is defined as the ratio of the actual pressure-recovery coefficient 
to the ideal pressure-recovery coefficient. Corresponding to Equations (8) and (9), there are 
two definitions for pressure recovery coefficient as follows- 


(3.11) and 


^ Sm 

3.9 PROCEDURE FOR DIFFUSER ANALYSIS 

For the purpose of analysis diffuser is considered as four separate diffusing passages and 
two mixing passages as shown in Figure 3.5. These passages are as follows: 

Diffusing passages-^ l-T, r-2, 2-3 A, and 2-3B 
Two mixing passages — > 3 A-4A and 3B-4B 

The essential element of the diffuser analysis is the evaluation of the performance 
of a single diffusing passage, including the determination of both the performance 
parameters and the outlet conditions. Using one or the combination of the following three 
options listed below can carry out diffuser performance analysis for a single diffusing 
passage (the specific passages for which they can be used are indicated in the parenthesis): 

Option 1: Empirical-Data Method (1-1', T-2, 2-3A, 2-3B) 

Option 2: Stream Tube Method (1-1', r-2, 2-3 A, 2-3B) 

Option 3: Mixing-Equation Method (2-4 A, 2-4B) 

Out of the above three options mentioned, the stream tube method has been used to carry 
out the diffuser analysis between Stations 1-3A and 1-3B. Option 3, i.e., Mixing-Equation 
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Method is always used in mixing passages 3A-4A and 3B-4B, regardless of the option(s) 
used up to Station 3 A and 3B. 

In the difliiser analysis, the Empirical-Data Method is not being used here due to 
lack of sufficient data and its less accuracy as compared the Streamtube method. 
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3.10 STREAM TUBE METHOD 

This method involves a stream tube analysis to calculate development of the velocity 
profile within the diffuser, as well as an analysis of the wall boundary layer to provide the 
location of separation points, if any. The model is indicated schematically below. Here the 
subscripts U and L refer to the upper and lower wall of the passage respectively. 



Figure 3.5 Stream tubes in the stream tube method calculation 
3.10.1 Assumptions 

The following assumptions specific to the stream tube method are made: 

• The flow can be represented by a number of stream tubes in each of which the flow 
is uniform and isentropic. 

• The boundary layer is identical to that of a flat plate with the same pressure gradient 
that exists in the diffuser and its effect on the inviscid flow is to produce an 
effective displacement of the diffuser wall equal to the displacement thickness of 
the boundary layer. 

• The stream slope varies linearly from one wall to the other at any axial location. 

• Separation of boundary layer is assumed to occur when the shape factor (i.e. form 
factor) H exceeds a critical value, Hsep which is an input to the program. 

• Subsequent to boundary layer separation, no further static pressure rise occurs in the 
section of the diffuser being considered. 

3.10.2 Input Information Required 

The input information required for a stream tube analysis is: , „ • 

1) The geometry of the diffuser in the form of radial position of each wall as function 

of axial distance. 

2) The following inlet conditions: 

i) Weighted mean total pressure, Pi. 

ii) Uniform stagnation temperature, Toi. 

iii) Mass flow rate, m . 
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iv) A velocity profile given by; 

-=%) 

U 

where U= an arbitrary normalization velocity 

f(y)= a tabulated function of against Y 

Y = nondimensional annulus height. 

v) The boundary layer displacement thickness, 6i . 

vi) The boundary layer shape factor, Hi. 

3) The number of stream tubes into which the flow is to be divided, N. 

4) An initial estimate of the boundary layer displacement thickness on each wall at 
every axial location of a geometric input point. 


3.10.3 Outline of Method of Solution 

The method of solution proceeds in the following steps: 

1) The determination of the static pressure at the inlet, and the mass flow, stagnation 
pressure, and velocity in each stream tube. 

2) The calculation of the sfatic pressure at each down stream location (specified by the 
location of the geometric input points), and the resulting radial location of, and 
velocities in each stream tube. The first time this calculation is performed, the estimated 
values of boundary layer blockages at all axial locations (supplied as input) are used; 
subsequent calculation are based on revised estimates of this blockage. 

3) The determination of the boundary layer displacement thickness on each wall, using the 
static pressure and wall stream tube velocities obtained in step 2. 

4) Repeat step 2 and 3 until the displacement thickness determined in step 3 are within a 
specified tolerance of those used in step 2. 

3.10.4 Inlet Conditions to the Stream-Tubes 

The inlet conditions for each of N stream tubes are determined as follows (the subscript j 

refers to a particular stream tube): 

1) Calculate the dimensionless coordinate Yj of the mid point of each (equal area) 
stream tube from the relations 



p 


J1 

^N+l 


1 


(3.13) 


where Rj+i = 


A. 


t+r/ 


% 
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2 ) 


3) 


4) 


5) 


Aj -R,") cos 

N ' ' ’ {^2 


R] = r + ■ 

l T_ 


COS 


t+Co' 


^Nh-1 “ 


V 2 
5 


u 


cos 


Cl+Cu^ 




At the compressor exit (Diffuser Station 1) it is assumed that 

1 = 1 


cos 


^C„+C: ' 


“i 


Find for each stream tube by interpolation from the velocity profile (tabulated 
as a function if Y) provided as input. 

Calculate pch, a characteristic density to provide a first guess for p from 


Pch “ 


R T„ 


-(3.14) 


where pi is the weight mean total pressure at the inlet. 

If the static pressure Pi at inlet is not given, make a first guess at this pressure from 


1 

Pi = Pi --Pch 


r • ^ 
m 




(3.15) 


where Ai = NAj 

Calculate the first guess at U from 


U: 


m 


Z^Poh Aj 


(3.16) 
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6 ) 


Using this value of U, calculate uj from the results of step 2 and the temperature Tj, 
from 


T. =T 

■^j -^01 



(3.17) 


7) calculate the density, pj, from 



A. 

RTj 


(3.18) 


8 ) 


calculate U from 

U = - 


m 




(3.19) 


9 ) 

10 ) 


If this value of U is significantly different from the previous value of U return to 
step 6. If not continue. 


Calculate the total pressure Pj in each stream tube from 


P; 


j _ 


^ T - '' 


v’^oiy 


T-I 


(3.20) 


11) Compute the mass flow in each stream tube from 

« 

m = PjUjAj (3.21) 

and the critical area of each stream tube Aj* from the usual compressible flow 
relations. 


A- JL 

a*“m 


r 2 ^^ 




1 + 

V 2 






y-l) 


(3.22) 


12) If the static pressure at inlet is calculated in step 4, calculate the new value of 
weight mean total pressure from 


(p,) =-ymV 

\ 1 /new • J J 

m j 


(3.23) 


where m, = mass flow in the jth stream tube. 
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13) If this new value of total pressure is significantly different from the given value 

of total pressure, calculate a new value of Pi from 

(^l)new ~ TTT (3 .24) 

VM/new 

and return to step 7. 

Hence, the quantities mj,Pj,Uj, Aj ,Pj are known at the inlet to each stream tube. The inlet 
dynamic head is calculated from 





j m 




(3.25) 


3.10.5 Determination of Stream-tube Properties at Downstream Axial Locations 

The conditions in each stream tube at each downstream axial location are determined in an 
iterative manner by estimating the static pressure at the axial location, computing the 
resulting total area of the stream tubes, and repeating the procedure until the computed total 
flow area equals the actual flow area normal to the passage center-line. The actual flow 
area is determined from the effective wall geometry obtained from the actual wall position 
and the estimate of boundary-layer displacement thicknesses. 

(a) Calculation Procedure; 

Subscripts, i refers to axial location of a geometric input point, and 
refers to a stream tube. 

i) As a first estimate, assume that the relative mid-stream tube streamline slope at i : 

P.26) 


where C - r 

bq bj 2 

= mid-stream tube streamline slope. 

For the first axial location downstream of the inlet, assume 



+ ^(CrU~CL) 


i+1 


since the relative streamline slope at the inlet is always assumed to be zero. 


ii) Estimate the pressure at the next axial location (i+1) from 
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1 



(3.27) 


Pi+i— Pi+Ap 

where TA, = effective total area at axial location i 


TAi 




2 

f \ 

2" 

- 


r„-5., 


r _1_ g 






^ " c +c 

cos-^^= — ^ 


> cos 

2 

- 

V 2 y 


^ 2 J 




iii) Calculate the static to total pressure ratio for each stream tube, Pi+i/P; and 
determine Tj/Toi > Oj , and Aj at i+ 1 each jfrom the usual compressible relations. 


iv) Compute total flow area based on initial estimate of streamline slope from 


(TaL=[|;Aj cosC,' ■ 

V j"' Ji+l 


(3.28) 


v) Re-estimate the relative slope of each mid-stream tube streamline from 




(3.29) 


ZK cosC,k)-^ cosC^ 
where Yj = 


Z(Aic cosC*) 


fc=i 


and k is a dummy index. 


vi) Compute revised total flow area based on new estimate of streamline slope from 
Equation (3.29). Repeat Step 5 until two consecutive values of (TA)'i+i are equal. 

vii) If(TA)'i+, does not equal TAi+i obtained from Equation (3.27), then repeat entire 
process from Step 2 with Pi , pj;, Aji,Cji in Equation (3.26) being replaced by the 

current values of pi+i, Pji+i, Aj ;^, , Cj.i+i • 
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viii) Continue the preceding process until ar® obtained at the exit 

station of the diffuser. 


(b) Determination of Boundary Layer Displacement Thickness 

With the velocity, u,, and density, pi, in the bounding stream tubes (j=l and j=N) known, a 
new estimate of the boundary layer displacement thickness along each wall is made in the 
following way; 


i) The momentum thickness at each axial Station i is computed from the following 
relation, using the momentum integral equation for boundary layers; 

. 6 / 





u.. 1 


" pji ^ 

% 

0.0076 v-^ r 4 7/ , 





+ 

n 11 /6 


(3.30) 


where j=l or N referring to L and U walls, respectively. And 

1 


-X;) 


COS 


4 


+ 






V 




(3.31) 


ii) The shape factor, H, of the boundary layer along each wall is computed from the 
relationship Dussord [30]; 

AHj4..=H„+70(H,,,-1.05)^ (3.32) 

Hj is never permitted to be less than 1.1, since this relationship is not valid for 
hi ghl y accelerating flows, and since the characteristics of the boundary layer when 
subjected to highly favorable pressure gradient are of little interest in the present 
application. Hj is also not permitted to exceed 3.5 as this value exceeds any 
reasonable separation value. 

iii) For all axial location at which H does not exceed Hsep (critical value of H at which 
separation takes place) on either wall, calculate the resulting displacement 
thickness along each wall from; 

For axial location where H>Hsep on either wall, proceed to step v. 

iv) Since this value of displacement thickness tends to interact substantially with the 
bounding stream tube velocity, a new estimate of the displacement thickness is 
determined from the following relation; 
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5j.i=5,j+0.2-^iiAL 


(3.33) 


hyercalLktioT^S ° estimate suppUed as input to the boundary 

y . 5 . . and 6 ^.i are determined from the solution to the following 

equations as shown below [35]; 

fe + l) 6-V„=S,ll-vJ+Vr„ (S.. +5.,) 

(vh, +i) 5"„,+V/Hi 5-V = 6-^,1 i-v,J+,;,„^ (5..,, +5«) 


where. 


Sji 

\|/ = iL_ 240 

■ ru-fr 


(Hoj,- 1.05) e,, 24 


x'i-x'w 7 


y,-70 


(H„.i-1.05) A9.,_, 




N.i-1 5o,i_, 

r i-i ~r 

U>‘-* U-l 

if no separation occurs along either wall at any axial location, proceed to Step vii. 

If H>Hsep on either wall at any axial location, set H=3.5 on the separated wall for 
reference purposes, and calculate the displacement thickness along the separation 
wall from, for example: 

k “S Cl +5i cos ^a.), = j(r„ cos C, -Ch cos Ci.. ' 

V ‘ 7t 

where, (TA)sep= effective flow area at the separation point 


gos Cl -^N cos ^ruf - (r, cos +5, cos 

r*r\ci K (C rL +C fUV 


cos cos- 


Cu+Cl 
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This assumes that the static pressure and the flow area remain constant after 
separation. 

vi) In the separated flow region, estimate the new value of displacement thickness 
from; 



vii) If the newly estimated values of Sj j are not equal to the originally estimated ones, 

^oj.i > return to the stream tube calculation with the newly estimated value. Repeat 
both stream tube and boundary- layer calculations until agreement is obtained. 

3.10.6 Determination of Diffuser Performance Parameters and Outlet Conditions 


The diffuser performance parameters are calculated in the following way: 


> Ideal pressure-recovery coefficient 

fA Y 

Cpi=l- 

UJ 


(3.34) 


> Actual pressure recovery coefficient 

P2-P1 

Cp 

qi 


(3.35) 


where qj is the inlet dynamic head. 


> Diffuser effectiveness 

(3.36) 

^p> 

> The blockage due to displacement thickness 


Eh=i- 


^U2 


"^U2 




^V2 


^L2 


(3.37) 


^ The outlet conditions of p2, p 2? n2(r) are directly available from the procedure, as well 
as p j , Uj , mj, Aj , for each stream tube. 
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> Pressure Loss for a Dump Diffuser 

The pressure loss for a dump diffuser was can be estimated in the range 
0.2<Ai/Ae<l by the equation(Barclay, 1972), 


^ = EXpi 



1 - 



+ 


r A ^ 
1 -^ 
A 




(3.38) 



Figure 3.6 Model dump diffuser 

For example, given the following dump diffuser data; 

Ae/Ai = 5.0, Mi =0.5, and/ = 1.4 . Using above equation (3.41) gives 
Pte/P,i = 0.854 


3.1 1 LIMITATIONS OF THE METHOD 

The major limitations of the method are as follows: 

i) Diffusers which are highly curved or in which streamline curvature effects are 
otherwise important cannot be treated. In these cases, the static pressure gradients 
across the passage are significant. 

ii) The results of the method tend to be quite sensitive to the inlet boundary layer 
properties; an accurate estimate of these properties at the diffuser inlet is difficult in 
practical situations. 

iii) The determination of location of separation is subject to error. 

iv) Dump diffuser analysis is difficult to carryout in this method as flow becomes 
turbulent. Assumption that the nature of flow being one-dimensional is no more 
valid in this region. Vortex formation and eddy losses are the dominant phenomena 
in this case. Two-dimensional analysis or the use of empirical correlation is capable 

of predicting dump losses. • • i j + 

The method does, however, represent a substantial improvement over the empincal data 

method. 
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3.12 MIXING EQUATION METHOD 

The mixing equation method can only be applied to the diffusing passages 2-4 A and 2-4B. 

3.12.1 Assumptions 

The assumptions employed in the method are: 

i) The losses in the diffusing process can be represented as the some of mixing losses, 
losses due to curvature of the passage, and expansion or contraction losses at the 
snout or cowl lip. 

ii) The mixing losses can be evaluated from the determination of the pressure recovery 
coefficient for incompressible flow. 

3.12.2 Input Information Required 

The information required in the mixing- equation method is: 

i) The geometry of the passage, in terms of inlet area, exit area, the total angle 
difference between the slope of the passage center-line at inlet and that at exit, and 
the mean radius of curvature of the passage center-line. 

ii) The inlet conditions defined bym, pi, qi, and Pi; the subscript 1 here refers to 
condition determined at diffuser Station 2. 

iii) The ratio of the area at diffuser station 2 occupied by the mass flow in the amulus 
under consideration to the annulus area normal to the engine center line at diffuser 
station 2', Ath/A 2 , i.e., for inner annulus 

(A,h) 2 A= tc [(r* - (rjA f ] 


where, (rtii) 2 A is determined from 

, ('ihhA 

m= ^Inrpudr 


where, m= mass flow in inner annulus 

r 2 A= radius of inner casing at station 2 A 
p, u = properties form diffuser solution at Station 2A 
A similar calculation may be performed for the outer annulus. 


3.12.3 Development of Mixing Equation for Incompressible flow 


(3.40) 


(a) Solution Procedure 

The method of solution proceeds in the following steps: 

i) Calculate diffuser performance between Stations 1 and 2, and the outlet conditions 
at 2, by either the empirical data or stream tube method. 
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ii) Detei^ne inlet conditions for the second stage of diffusion in a form suitable for 
the diffuser analysis methods to be used between Stations 2-3A and 2-3B and the 
inlet conditions to the snout (if any). 

iii) Determine the stagnation pressure on the dome or cowl. 

iv) If the diffuser has no snout, determine the inlet conditions to the annuli by assuming 
an isentropic expansion around the dome. 

v) If the diffuser has a snout, calculate the diffuser performance between Stations 2-4A 
and 2-4B and the outlet conditions at 4A and 4B, by using the empirical data 
method or stream tube method between Stations 2-3A and 2-3B and a mixing 
analysis between stations 3A-4A and 3B-4B, or by using the mixing equation 
method between Stations 2-4 A and 2-4B. 

When the stream tube or empirical data methods are used between stations 2-4, it is always 
implied that the methods are used between Stations 2-3 followed by a mixing analysis 
between Stations 3-4. 

For diffusers without snouts (or cowl), the analysis of the flow between Stations 2-4 
is independent of the option employed between Stations 1-2. 

(b) Inlet Conditions for the Second Stage of Diffusion 

For diffusers without snout or cowl, the only additional information required is the mass 

• • • 

flow split, mA, ms, ms. This is supplied as input from the airflow subprogram. 

For diffusers with snout, the information required in addition to the mass flow split, 
consist of the flow properties Ah, Ua, p (or E) to be associated with the flows into the inner 
annulus, the outer annulus, and the snout or cowl. Using stream tube method these can 
be obtained, between Stations 1-2, as below: 

i) The area of the passage at station 2, Ah, 2 A and Ath, 2 B , which are occupied by the 

« • 

mass flow mA and ms are determined from Equations 

I . (.nhhA 

(Ati,) 2 A= 7ti(ra,)2A^ -(r 2 An \2nrpudr (3.41) 

applied to both mass flows. Then Ath^s is determined from the equation; 

Aih, 2S = Aa -Ath, 2A -Ah, 2 b (3 .42) 

where, A 2 is the total mass flow at the combustor inlet. 

ii) The inlet area average velocities are determined from 

u,„=^judA (3.43) 

^th,2A 0 

and a similar equation for Ua, 2 B . 

iii) The inlet velocity profile parameters are calculated from 
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-^lh. 2 A 

pdA 

(32a= — ^ — (3 44) 

'^th,2A '^a.2A 

and a similar equation for ^2b. It is assumed that p2s=l. 

For the purpose of the difluser analysis between Stations 2 and 3, the areas of 
diffusing passages at Station 2' are assumed to be 


Ath.2A~ A<h, 2 ACOsqc,A (3.45) 

and 

Ath, 2 B = Ath, 2 B COS qc,B A (3 .46) 


where qc angle of the slope of the passage centerline relative to the axis of the 
combustor. This assumes in effect that the inlet areas to the annular diffusing 
passages are the areas occupied by respective flows at Station 2, projected normal 
to the passage centerline. 

iv) The inlet dynamic head at the diffusing passages are calculated from: 


1 m^jJ^u^ 


12A 


^th,2A 


(3.47) 


and a similar equation for q^g . The inlet dynamic head at the snout or cowl is 
determined from; 


q 2 s = 


1 m^ Uro2 
^ Allies 


( 3 . 48 ) 


(c) Conditions on the Dome 

If the diffuser has no snout or cowl, the stagnation pressure on the dome is assumed to be 
constant and is determined by 

P.=P,+(l-K,)i„ (3.49) 

where 92 =rP 2 ( 9 mi) 

The coefficient K4 is an input quantity representing the number of dynamic heads 
lost in the mixing and flow processes occurring around the dome. The magnitude of Ka will 
depend upon the shape of the dome and the effectiveness of the diffuser between Station 1 
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and 2, and is difficult to determine a priori. It is expected that its value would be in the 
range 0.2 - 0.5. 

For diffusers with snout (or cowl), the static pressure just inside the snout, P 2 S', is 
determined by assuming a sudden contraction or expansion process between the areas Ath, 2 s 
and A 2 S' (the actual geometric areas of the snout entrance) 


*I2S' ~ 


K '' 

^th,2s 

V ■^2s' J 


I2s 


(3.50) 


The total pressure on the dome is then calculated from 


Pd=P2s'+(l-Kd)q2,. A (3.51) 

where Ka is an input quantity representing the number of dynamic heads lost. 


(d) Diffuser Performance between diffuser outlet and Annuli diffusion region 
(Stations 2 and 4) 


The procedure for determining the performance of the diffusing passage between Stations 2 
and 4 depends upon the diffuser analysis method employed. Using stream tube method 
the procedures are as follows: 

In the stream tube method, three steps are involved since the basic stream tube 
analysis is only applicable between Stations 2 and 3, and do not consider any loss due to 
expansion or contraction of the flow at the snout lip. 

The steps are as follows; 

1) A basic stream tube analysis, as previously described, is performed between Stations 

2A-3A and 2B-3B. The inputs required for this analysis are m, Aih, 2 A', P 2 , U 2 (r), , 

H 2 A and similar set for the outer annulus. In addition it is assumed that the boundary 
layer properties on the snout lip are 5=0 and H=1.4. The only calculated quantities 
which are used subsequently are the static pressure, Psa and Psb and the profile 
parameters, Psa and Psb which are computed from Equation (3.47) applied to the 
Stations 3 A and 3B. 

2) It is assumed that the profile mixes to a uniform one between Stations 3 and 4. The 
static pressure at Station 4 in the absence of area mismatch at the snout lip is 
accordingly obtained from Equation: 


P' -P 

^ 4A _ . 


1 


I+A3A 




HAy 


P 


3A 


••3A 


‘•4Ay 


(3.52) 


fisA P3A 

3) The effects of any area mismatch at the snout are cmdely approximated by assuming 
that the total-pressure loss between Stations 2 and 4 due to this process is given by the 
sudden expansion or contraction relations, and that this loss may be applied directly at 
Station 4; i.e., for example, in the inner annulus: 
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P -P' 

4A 4A 
^2A 


1- 


( k ^ 


■^th,2A' j 

'v ^2A' > 


(3.53) 


if 


*^2A' 


‘^th,2A' 


< 1 , then P' 4 A is the total pressure corresponding to P'aa, To, , and Aaa. 


3.13 LIMITATIONS OF THE DIFFUSER SUBPROGRAM 

The analytical methods used in the difiuser subprogram contain several limitations; 

^ Flow in the diffuser is assume to be unaffected by downstream conditions. If the 
flow is well matched at the cowl, this assumption is expected to be valid. 
Otherwise, serious errors will undoubtedly result. Cowl design should, therefore, 
receive careful attention. 

> No provision is made for including splitter vanes in the diffuser analysis. 

> The program does not attempt to treat flow after separation. 

3.14 STRUCTURE OF THE SUBPROGRAM 

The diffuser subprogram consists of 

> A subroutine which sets up starting conditions and directs the diffuser calculation 
through the other subroutines in the appropriate sequence. 

> A group of subroutines which perform the stream tube analysis. 

> A group of library subroutines which provide gas properties and interpolation 
procedures. 

Subroutines and flowchart of the diffuser subprogram is given in the Appendix A2. 



CHAPTER 4 

AIR-FLOW ANALYSIS OF MAIN COMBUSTOR 


4.1 COMBUSTION CHAMBER FLOW PATH ANALYSIS 

Analysis of the combustion chamber flow path is closely related to, and proceeds in parallel 
with, the inlet diffuser analysis. These analyses must be coordinated because the combustor 
cowl and passage contours are very important to successful diffuser operation. Conversely, 
diffuser pressure recoveries must be known in order to select appropriate cooling and 
dilution hole sizes to get the specified flow distribution. 

For most applications, the preferred flow path is the shortest one that meets all design 
requirements. Increased length adds weight to the engine and requires more liner cooling 
flow, which reduces available combustor air. 

4.2 PERFORMANCE PARAMETERS 

The main performance parameters for an annular combustion chamber are (see Figure 4. 1): 

^ Combustor dome height (Ha): the dome height of the combustor is important 
because it has a significant effect on altitude relight capability. The dome height 
must be large enough to permit the combustion system to satisfy the altitude relight 
requirements of the engine. 

^ Combustor dome velocity (Vj): the combustor dome velocity is defined as the 
average velocity for all of the combustor dome flow, immediately downstream of 
the dome, between the inner and outer combustor liners. In general, combustion 
flame stability is reduced as dome velocity is increased. 

> Combustor length to dome height ratio (Lc/Ha): all combustors must have 
sufficient volume and length to accommodate a low-velocity flame stabilization 
region and a higher velocity mixing region, where the hot combustion products are 
mixed with the excess dilution air. The total necessary volume depends on space 
rate, residence time and other loading parameters, but the overall flow pattern 
required is well established. The desired flow patterns determine, to a large extent, 
the overall shape of the combustor, as defined by the combustor length to dome 
height ratio (Lc/H<j). Since the basic flow patterns can be scaled to larger and 
smaller sizes, Lc/Hd is essentially independent of engine size. Therefore, Lc/Hd 
provides a useful means of comparing the combustors of different engine models. 

^ Combustor passage velocity (Vp): liner passage velocity for a combustion system is 
the velocity of the flow between the inner and outer annulus of the system. The 
passage velocities depend on the mass flow rate and passage flow area, but it is 
defined as the velocity at the dome exit plane with the total passage air flow. The 
passage velocities should be low enough to provide a uniform flow to the 
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combustor with high static pressures and low total pressure losses in the passages, 
and high enough to provide good convective cooling of the combustor liner. 




Figure 4.1 Combustor chamber variables 

^ Fuel injector spacing (B): ideally, for good combustion system performance, the 
fuel should be injected at an optimum number of points. However, practical 
considerations usually limit the number of fuel injectors to a few only than optimum 
number. Fuel injectors along with the swirl cups which are used with each injector, 
are costly and heavy. Also, a large number of injectors increases maintenance time 
and costs. The ratio of combustor length, Lc, to circumferential spacing of fuel 
injectors, B, is a critical design parameter. If the fuel injectors are too widely spaced 
and the Lc/B parameter is too low, high-temperature regions may appear 
downstream of each fuel injector at the turbine inlet plane which may cause the 
pattern factor to be too high to meet requirements. 

> Space rate (SR): the combustion system space rate is a measure of the fuel energy 
released per unit volume of the combustion chamber. 

> Reference area (Rref): it is the area at the point of maximum flow area between the 
combustor casings. 

> Reference velocity (V^f): it is a measure of an “average” velocity through the 
entire cross-section area between the inner and outer casing walls of the combustion 
chamber. This velocity affects the residence time of fuel-air mixtures in the burner 
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and the basic flame stability of the system. Although the reference velocity does not 
describe a physical velocity at any particular place in the system, this parameter 
provides ^ a convenient method of comparing different designs. The reference 
velocity is calculated based on compressor exit flow, and air density calculated 
using compressor exit total pressure and total temperature. 

Reference velocity can also be defined in terms of the Tnaxirmim combustor 
cross-sectional area. 

^ Reference velocity head (<iref): the reference velocity head is the difference 
between total and static pressure at the inlet air density and velocity. 


qref- 'A(P3 V^tef) 


The magnitude of Qref relative to pressure drop across the combustor dome 
and liners is important to uniform air distribution and good dilution jet penetration. 

> Inlet velocity head (qj): it is the difference between total and static pressure at the 
compressor discharge. 

firef = 'A(P3 V^a) = Pt 3 - Ps3 

4.3 AIRFLOW DISTRIBUTION 

Combustor airflow is distributed to different combustor locations to achieve different 

design goals. Flows are generally specified as a percentage of total combustor airflow. 

The airflow distribution does not vary significantly with combustor operating conditions. 
The most important airflows are; 

1) Compressor exit flow. 

2) Turbine cooling airflow which bypasses the combustor. 

3) Combustor air flow (1-2) 

4) Fuel atomizing airflow admitted through the fuel injectors to the fuel into small 
drops. 

5) Swirler airflow admitted through a swirler around the atomizer to provide a strong, 
well-mixed recirculation zone within the primary zone. 

6) Primary air jets which interact with the swirler flow to close the primary, zone 
provide rapid fuel-air mixing and oxygen to complete combustion reaction in the 
secondary zone. 

7) Dilution air jets down stream of the primary zone to provide the desired exit 
temperature profile. 

8) Dome cooling airflow. 

9) Liner cooling airflow. 

10) Combustor dome flow (4+5+8) 
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4.4 DOME AIR FLOW 

Dome flow consist of fuel atomizing flow, swirler flow and dome cooling air flow 

a. Fuel ^ atomizing flow, the atomizer airflow should be at least two to three times 
atomizer fuel flow. Lower levels of atomizing flow will reduce the atomization. 
Higher levels will tend to reduce lean blow out capability. 

b. Swirler flow: the swirler is used to set up a redrculation zone for stability and to 
provide good fuel/air mixing for low smoke. For the sake of low power combustor 
efficiency and lean stability, a swirler flow level much higher than the nominal 11% 
is not recommended. 

c. Dome cooling flow: the last component of dome flow is dome cooling air flow. 
Dome cooling flow is estimated based on the total area to be cooled (dome area less 
exit area), and a cooling rate based on previous experience. For a typical combustor, 
dome-cooling flow will be 10 - 15% of combustor airflow. 

4.5 COMBUSTOR LINER AIRFLOWS 

The combustor liner flows are mainly through primary, secondary and dilution holes. The 
primary holes are typically located at an axial location about one dome height downstream 
of the dome. Primary air serves two purposes. First, part of the primary air interacts with 
the atomizing and swirler airflows and a portion of the dome cooling air to prevent smoke 
formation in the primary zone. Second, all of the primary air, along with all atomizing, 
swirler and dome cooling air is needed to reduce the overall equivalence ratio down stream 
of the primary holes to less than stoichiometric, thereby providing oxygen to complete the 
combustion reactions initiated in the primary zone. The zone immediately downstream of 
the primary zone is called the secondary zone. Secondary hole airflows ensures the 
complete combustion of the combustible mixture that was burnt in the primary zone. The 
remainder of combustor airflow is available for liner dilution holes down stream of the 
secondary zone. Dilution airflow helps in getting desired temperature profile at the 
combustor outlet or turbine inlet. 

4.6 DETAILED FLOW ANALYSIS 

As a minimum, a detailed one-dimensional flow analysis is needed to determine pressure 
distributions within the liner passages and combustion chamber. Total and static pressure 
distributions throughout the system are required to determine dilution and cooling hole 
sizes necessary to obtain desired airflow distribution. 

For one-dimensional analysis, combustor is divided into a large number of short 
segments. The start of each segment is defined by two axial positions or stations as shown 
in Figure 4.2. The first station is normally at combustor dome. An additional station is 
located on the centerline of each row of liner dilution or cooling holes. 
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Figure 4.2 Layout for one-dimensional combustor analysis 

The basic compressible equations are solved at each station. The initial pressures 
are based on passage recoveries calculated on diffuser analysis. Between calculation 
stations, the effects of friction losses and momentum exchange between the dilution jets 
and internal flows are calculated and results are reflected in calculations at the next station. 
An interactive solution procedure is used to balance flow and pressure distribution. 

4.7 ANALYTICAL METHODS AND PROGRAM DEVELOPMENT 

4.7.1 Objective of the Airflow Subprogram 

The main objectives of the airflow subprogram are: 

> To calculate the air mass flow and pressure distribution in the flame tube and 
annuli and, hence, the overall pressure loss. 

> To provide a starting point for the heat transfer calculation by working out the 
temperature and velocity distribution of the gases in the flame tube and annuli. 

To facilitate analysis, flow in the combustor is broken down into distinct but interacting 
streams; 

1 . The annulus air (inner and outer annuli). 

2. The air flowing into the flame tube through: 

(a) Swirler. 

(b) Dome holes. 

(c) Penetration holes in the flame tube. 

(d) Wall cooling slots in the flame tube. 

3. The mixture of fuel, air, and combustion products flowing through the primary 
zone (the flame tube up stream of the secondary holes). 

4. Parallel streams flowing through the flame tube down stream of the primary zone; 

(a) The hot stream, a mixture of combustion products, unburned fuels, and air. 
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(b) The cold streams, those portions of the penetration jets and wall cooling 
iilms, onginating upstream, that remain unmixed at a given axial position 
The following sections give a detailed description of the calculation methods. 

4.7.2 Assumptions 

The air flow calculation rests on the following main assumptions; 

• The primary zone can be treated as a stirred reactor: 

a) Static pressure, temperature, and fuel air ratio are uniform through out the zone. 

b) Mixing and burning occur instantaneously. 

This is one of the major simplifications in this work. It is particularly invalid at low 
pressures. 

• The flow can be considered to be one-dimensional. This represents a considerable 
simplification, particularly in the case of the gases in the flame tube. 

• The rate of fuel burning as a function of axial length can be specified as inp ut This 
assumption makes the calculations easier but a realistic estimate of the fuel burning 
rate will be difficult to make, particularly (as mentioned above) at low pressures. If 
the fuel air ratio at any ^ial position, which will not be known in advance, exceeds 
the stoichiometric value, the excess fuel is assumed to be available for burning at the 
next downstream position. 

4.7.3 Calculation Procedure 

In calculating the flow conditions throughout the combustors, the air flow subprogram 
goes through the following routine for each iteration on the mass flow split: 

1 . The mass flow rate through the dome is obtained from the flow split, and the total 
pressure on the dome is obtained from the diffuser calculation. 

2. The pressure drop across the dome to the primary zone is found by equating the 
combined swirler and the dome hole flows (functions of pressure drop) to the total 
dome flow. 

3. With the primary zone static pressure established, flow conditions in the aimulus and 
the flow through flame tube holes into the primary zone are computed for calculation 
stations up to secondary holes (the end of primary zone). 

4. The portion of the secondary hole flow that recirculates into the primary zone is 
calculated, either directly by the program or from a fraction specified as input. 

All of the contributions to the primary zone air flow (through the swirler, dome holes, 
and flame tube holes up to the secondary holes) are now knovra. In addition, the fuel flow 
into the primary zone is known from input. 

5. The temperature rise due to combustion is calculated for the primary zone treated as a 
stirred reactor. (Excess fuel above stoichiometric is carried downstream). Flow 
conditions at the end of the primary zone are then found. 

For the remainder of the calculation the combustor is divided into axial control volumes 
bounded by adjacent calculation stations k and k+1. Calculations proceed down stream 
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from station to station. In general, the flow conditions at the upstream station, k, are 
known from the results of the previous calculation step. Conditions at the downstream 
station are found as follows: 

6. From the ^nulus to flame tube pressure drop at station k the mass flow rate through 
the holes just downstream of k is found. The mass flow re maining in the annulus is 
then know. 

7. For the annulus, the equations of continuity, momentum, energy, and state are solved 
for the pressure, temperature, velocity, and density at station k+1. 

8. Correlations for flow and mixing of wall and penetration jets are used to compute the 
fluxes of mass, momentum and enthalpy contained in the jets entering the control 
volume (through holes at stations k and in the form of residual jets originating 
upstream of k) and leaving the control volume (residual jets at k+1). The net transfer 
of mass, momentum, and enthalpy from the jets into the main gas stream in the flame 
tube may thus be found. 

9. The equations of momentum, energy, continuity, and state are written for the hot 
stream in the flame tube (excluding residual jets) between stations k and k+1. Heat 
addition is included for fuel burned, up to stoichiometric; any excess fuel is carried 
downstream. The solution of these equations yields the temperature, pressure, density, 
and mass flow at station k+1. 

The sections that follow describe the detailed calculations that are performed in each of 
the above steps. 

(a) Flow through Holes and Swirler 

The mass flow rates through the holes in the dome and in the flame tube walls are 
calculated from the discharge equation; 

=CaAh(2p„(p,„-Pft))" (4.1) 

where Pan = total pressure in annulus 

”■ Pan qan 

qan~ PanU an/2 

Cd = discharge coefficient 

= 0.6 for dome holes; for other holes types the data are draw from 
[10], [11], [12], [13], [14], and [15]. 

(b) Pressure- Ratio Correction 

The effect of flame tube velocity on discharge coefficient can be predicted by [14]; 

Cd = Cd^ (0.75+0.25 pan/pft _ — _(4-2) 

where Cd = actual discharge coefficient which is a function of pressure ratio. 

Cd^= corrected discharge coefficient which is independent of pressure ratio. 

(c) Swirler Flow 

The equation used in the program to determine the flow through a swirler is given [16]; 
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where 


AP 

sw _ 

Qtef 



sec^ P 



m! 


m: 


ref 


(4.3) 


A = area 
m = mass flow 
q = dynamic head 
P = swirler angle 
sw = swirler 

ref = reference conditions 
ft = flame tube 
Ksw = constant 

— 1.3 for straight bladed swirlers; and 1.15 for curved wall swirlers. 


(d) Recirculation 


In an actual combustor a complex pattern of recirculating flow exists, particularly in the 
primary zone. The fraction of secondary hole air flowing into the primary zone, msy,p/msy, 
may be specified as input to the computer program; msy,p/msy = 0 if not specified as input 
to the program. It may be generated in the program from the following empirical formula 
[17]: 

^ = 0.5 sin 

tttgy 

The air entering the primary zone by recirculation from the secondary holes is treated in 
exactly the same way as air entering through swirler, dome and wall holes. Air going 
down stream is treated as ordinary penetration jet air. 


it/z 


(4.4) 


(e) Primary Zone of Flame Tube 


The pressure in the primary zone, pi, is found from the pressure in the diffuser upstream 
of the dome, obtained by the diffuser subprogram, and the pressure drop across the dome, 
which is a function of the flow through the dome and the characteristics of the swirler and 
the dome holes: 


m. 


■Pi =• 


2p 


ref 


■^dh^dh ( 



sec^ P,, 


where ma = total flow through dome 
■“ msw mdh 

mdh = flow through dome holes 

Cdh = discharge coefficient of dome holes. 

Ksw= coefficient in swirler pressure drop equation. 

Psw = blade angle of swirler 

Aft= cross-sectional area of primary zone 


(4.5) 
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(f) Temperature Rise 

The energy equation for combustion in the primary zone is given by (neglecting the 
velocity terms as mean velocities in the primary zone are quite small); 

+ q = (m, + mjb) (4.6) 

where mi = air mass flow rate 
mfb = fuel burning rate 

q = heat release rate 
hair= enthalpy of air 
hprod = enthalpy of combustion products 

The subscripts 1 and 2 refer to conditions before and after combustion, respectively. 

The enthalpies of air and combustion products are calculated from the expression 
(neglecting the dissociation effects): 

h = |cpdT (4.7) 

where the specific heat of air as a function of temperature is correlated by: 

Cp =0.2419-0.8181xl0'^T+17.91xl0’®T^-2.743xl0'’^T^ (4.8) 


The specific heat of combustion products is given by [18]: 


Cp =0.2419 + .103f-(0.8181-22.6Qxl0''T + (17.91-29.6f)xl0'®T^-(2.743-.35f)xl0-''T' 

(4.9) 

where f = local fuel air ratio. 

This equation is suitable up to stoichiometric values of the fuel-air ratio which can be 
obtained from the relation; 


stoich 


0.0867 


1 + H/C 


J-i-3H/C^ 

where H/C = fuel hydrocarbon ratio by mass. 


(4.10) 


The combustion of fuel in a vitiated gas, when the enthalpies don’t include the effects of 

dissociation, the heat release rate, q, is adjusted to account for dissociation and 
recombination. In general, for combustion between conditions 1 and 2; 


m^hp +3x10 


-26 


m^ T, ■ m^T; 

1 + f, 1 + f 2 


(4.11) 


where m^ = rate of burning of fuel 


f ^'hrw 

iilliL. 


ararrfc?! 
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For combustion in the primary zone the temperature, T 2 , is very large and the term 
m Equation (4. 1 1) may be ignored. Equation (4.6) is now solved for T 2 , the primary zone 
temperature. 


(g) Density and Velocity 


The density is found from the equation of state: 
P2 


P2 


RT, 


and the velocity from the continuity equation: 

m. 


U2 = 


P2A2 


(4.12) 

(4.13) 


where m 2 = mi 


(h) Hot Loss 


The “hot loss” of total pressure, which is the change in total pressure due to heat addition, 
is given by: 

AP = P, -Pj 



1-^ 


P ^ 

PW 


(4.14) 


4.7.4 Flow in Annulus 


The equation of continuity, momentum, energy, and state are set up for a section of one 
annulus between adjacent calculation stations, shown as 1 and 2 in the sketch below 
(Figure 4.3): 



Figure 4.3 Calculation stations in annulus 

The resulting equations are then used to calculate conditions at Station 2 from the 
conditions at Stationl. 
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(a) Continuity Equation: 


The continuity equation for flow in the annulus between Stations 1 and 2 is written as 
follows: 


Pan.l^ an.l'^an.l^ an,2 


m 


an,l 


Pan,2^an,2 “^3:1,2 


(4.15) 


(b) Momentum Equation 

Due to the presence of annulus wall holes between stations 1 and 2, the equation for the 
conservation of momentum becomes complicated. To deal with this complication the 
calculation is carried out in three steps. 


1) The flow through the hole, mn, is found from the hole type and size, the pressure 
drop across the hole, and discharge coefficient data. 

2) It is assumed that the air remaining the air remaining in the annulus at station 1', 
immediately downstream of hole, undergoes a sudden expansion. The resulting 
total pressure loss given by the following empirical expression [19], from which 
Pan,!' the pressure at station T may be calculated: 


Paa.-Pan,r =1-85 


Pan,l U 


2 r 


an,l 


m. 


1.36m, 


X5+.242M.,i“’) 


(4.16) 


With the total pressure at Station 1' known, the other flow conditions at Station V 
may be found from the equations for the isentropic flow of a perfect gas between 
Stations 1 and 1'. The approach is first to calculate the Mach number at Station 
r, Man.r from the following equation: 


^an,r _ ^an.2 


Aan,r ^ 


aikl 




M.,,. 


rJ 


•I Y“1x4- 2 


7+1 


(4.17) 


2(H) 


where To = stagnation temperature and assumed to be constant across the 
expansion 


u. 


T + 

an.1 2 C 


The static temperature and static pressure are then determined from. 


T . 

■*‘00,1 __ 




-1 
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Then from the equation of state the density at T is: 


Pan,l' 


Pan,!- 

RT, 


an,r 


From continuity equation, the velocity is: 


m 


Uan.!' = 


an,r 


Pan.1- A 


an,r 


all of the 


Thus, the assumption of a sudden expansion at Station T permits 
necessary flow quantities at that station to be completed. 

3) Now apply the momentum equation for flow between Stations 
r and 2; 

A„, u.,, f p. ^dx -' f u 

where F= fanning friction factor PO] 

= 0.0035+0.264 (Re)a„,i 
Aw = Wetted wall area per unit length. 

A simplified form of equation (4.18), on ignoring small pressure losses due wall 
slope and friction, is as follows: 


an,2 

(4.18) 


1 ^^^P"U - g !dx = + Cb)p^,.U3^,"(x, - X,) 


■(A19) 

^^ 2 2 

where, c'a and Cb are the inner and outer wall areas/ unit length between 1 and 2. 

The integral pressure term can be evaluated if it is assumed 

Stations r and 2 is incompressible and that the mixing rate ^ ^ J 

diffiision. With these assumptions the integral pressure term may be integrated as 

follows; 

The integral pressure term in the annulus momentum equation (Eq. 4.18) may be 
integrated by parts to get; 
dA 


•> dX ^^m,2 
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Now, the final momentum equation thus becomes; 


Aan,l + Aan.2 



-p„2)-f(c» 



-u 


an,!' 


) 


(4.20) 


(c) Energy Equation 

The energy equation for flow in the annulus between stations 1' and 2 is; 




man,2 + 


f 

JqwCAdx = 




h(T„,)+^ 


m 


an,2 


where h(T)= static enthalpy of annulus air at temperature T 

q^^ = heat transfer rate fi'om unit area of the flame tube wall 
to the annulus air. 


The casing temperature is normally assumed to be equal to the compressor delivery 
temperature because the heat transfer rate between casing and annulus air is assumed to 
be negligibly small. 


Since the variations in the annulus air temperature are normally small, h (T) may be 
assumed equal to Cp (T-Tb), where Cp is a mean specific heat. If mean values are assumed 


m 

for and Ca, the simplified energy equation becomes: 


u 




an,l 


CpTjn,] + ^ 


/■ 


m^2+qwCA(x2-x,) = 




2^ 


^an,2 ^ 2 


m 


aii,2 


(4.21) 


(d) Equation of State: 

The equation of state for the annulus air is; 

Pan.1 _ Pan,2 _ ^ (4.22) 

Pan,l^an,l Pan,2^an,2 

4.8 METHOD OF SOLUTION OF ANNULUS EQUATIONS 

Since the mass flow through the holes between stations 1 and 2, mn, is known, the down 
stream flow in the annulus, man, 2 is simply. 

man.2 

There remaining four unknowns, Uan .2 , Pai .,2 , Pan ,2 and Tan ,2 , which can be found by 
solving equations (4.15), (4.20), (4.21) and (4.22). 
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These equations can be easily reduced to a quadratic in Uan 2 : 


B 2 V2'+B, U3^,-B„=0 

c„ A.„,+A., 


(4.23) 


where B, =1 — 


R A 


an,l 


B, =2 
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r. ™,2 , 

Pan,!' V 
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*an,r 


^an,l + 2 C 


an,2 


Tan,l + 


^an.1 ^ 

V ^^an.1 J 


1 ^ (^A ^3X^2 ~ Xi-) 


2A. 


an,! 


q„ C^(x,-x,) 


P an,2 


The solution of Equation (4.23) gives, 

„ _-B,+VB|^+4B„B, 


2B, 


(+ve root) 


(4.24) 


With Uan ,2 known, Equations (4.15), (4.20), and (4.22) are solved for , p^i , and 

Taii,2 • 
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4.9 CHARACTERIZATION OF JET FLOW AND MIXING 

The flow characteristics of jets entering the flame tube through holes and cooling slots or 
rings is treated by the use of empirical correlations for: 

> The length of the jet (measured along its center line). 

> The cross-sectional area of the spreading jet as a function of jet length. 

> The jet center-line velocity and temperature as functions of jet length. 

> The development of the non-dimensional velocity and temperature profiles across 
the jet with jet length. 

Two sets of correlations are provided, one for wall jets (where the air enters the flame 
tube through a cooling slot or ring in a direction parallel to the wall) and one for 
penetration jets. These two types of jet are shown below ( Figure 4.4); 
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Figure 4.4 Combustor penetration and wall jet. 

The characteristics of penetration and wall jets and of the three mixing models (any one 
of them can be used in the program) are as follows: 

4.9.1 Penetration Jets: 


The dimensions used in the describing jet behavior are shown in the sketch below: 



Figure 4.5 Penetration jet parameters 

Where x, is the direction of flow of main-stream gas in the flame tube 

y, is the jet penetration distance normal to the flame- 
tube wall 

s, is the distance along the jet center line 

T], is transverse width of the jet (measured normal to the jet center line 
along lines of constant s), and 
(p, is the angle between the initial jet trajectory and the 
flame tube wall at the entry hole. 

(a) Initial Jet Angle Data 

The jet mixing correlations, and the momentum equations in the flame tube, make use of 
the initial jet anglecp. This is available as a function of pressure-loss factor for different 
hole types in the library data. Where data do not exist, the following correlation is used: 
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sin> = -^ (4.25) 

^dn 

where Can is the asymptotic value of Ca as 1 + tends to infinity. 

^an 

This approximation gives reasonable results for round flush holes. 

(b) Penetration Distance: 

The penetration distance y, measured fi'om the wall to the jet center line, is calculated by 
an empirical correlation [21] through [23]. 


y = 0.87d, 



0.47 




(4.26) 


where dh = effective axial length of hole 
=(Ca)’^^ X axial length 
Uj,o= initial jet velocity 

p^C<jXhole area 

It is assumed that there is no change in density of the jet as the air passes through the 
hole. 
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(d) Transverse Velocity Profile; 


The correlating eq^tion for the jet velocity profile, expressed in terms of the transverse 
coordinate, T|, are obtained from [22]: 


Uj,,=Uft + 


■ 

( \ 

- 

0.5 cos 

71 ri 

+ 0.5 


^ ^ ^1/2 ^ 



(u j.cL ~ tift ) , for 0 < •q < 


u.,,=u,+1.74e ^(uj.cL-Uft) , 
where, from [22]; rii /2 =0.215s' 


for T]>ri 


Yl 


(4.30) 


(e) Transverse Temperature Profile: 

The transverse temperature profile is assumed to be similar in form to the transverse 
velocity profile. The expressions relating temperature to velocity are: 


'^j.CL “ Tft u - Uft 


(4.31) 


'^j.CL.s '^it _ 4.8 


T - T 

^■j.o •'■ft 


6.5 


^j.CL.s' ^ft 
V ’^i.o-Uft J 


(4.32) 


where 'n"=Ti/1.4 

The constants are taken from Forstall and Shapiro [24], and Squire [25] 


(f) Calculation of Jet Properties: 

The cross-sectional area of a jet in a plane normal to its center line may be represented by 
[24]: 


A,,„=0.12q*" 

where ti* = the value of tj at the assumed jet boundary. 

The jet area normal to the combustor axis is then given by: 

0.12 q*' ® 


A.=- 


where 


COS<i3j 

(pj = local jet angle 


COS 


1 + 




^dy 

Vdxy 


- 1/2 


(4.33) 


(4.34) 


(4.35) 
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Other jet properties needed for solution of the flame tube equations may 
follows: 

Mass flow rate, mj = JpjU.^O.Uxl.ST)''* dt] 

0 

Axial component of jet momentum flux, 

T|" 

0 . 12 X 1 . 811 *'* dTi cos^. 

0 

Enthalpy, h- = JpjUj., 

0 


'=,(Ti,-Tj+in.; 


0.12xL8Ti°'dTi 


be calculated as 


(4.36) 


-(4.37)) 

~(4.38)) 


4.9.2 Wall Jets: 

The symbols used are defined by the Figure 4.6: 


r " ^ 




s. 

'S. 

•v, 

; y' »■ » - 

“ft 


Figure 4.6 Wall jet parameters 
(a) Maximum Velocity: 

The maximum velocity in the wall jet is given by [26] 

When (i)x = 0, s = 0 


(ii) X = X + dx , 


Uj.CL=% 


s = s + 




U; 


dx 


i.CLj 


where 


^O.mjy, 

" v(l2s')- 

Uh = initial jet velocity 

Pan^d X slot area 


(Uh-Uft) 


(4.39) 


a = 


1.06 

1-Uft/Uh 



and 
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The above equations apply for values of Uj,cL between Uh and Uft . Otherwise, Uj,cL is set 
equal to Uh. 


(b) Velocity Profile: 


For analysis purposes, the velocity profile is divided into two sections about its 
maximum. 

i) Near the wall ('n<5i) the equation from [26] is used: 


j.n 

j.CL -Ufl 



il<5i 


(4.40) 


where 11 = perpendicular distance from wall 
C = 1.0 
n = l0 

5, = 0.0109s' 

Since n is so large, we can ignore the precise shape of the profile and to assume: 


dTi = 0.955, (4.41) 


r U ■ — Uft 


0Uj.CL-Uf 


ii) 


In the outer region (r|>5i) the equation from [27] is given by: 

Uj ,, - Uft _ .k ^'2 

Uj,CL-Uft 


(4.42) 


where k = 0.693 

, q-S, 

'n = — — — 

qi/2-5, 

= 0.065s' [36] 

(c) Calculation of Jet Properties: 

The equations for the properties of wall jets are similar to the penetration jets from 
Equations (4.43-4.38). The major difference is in the expression for elements of area. 


Aj 


mj 



( on the inner wall) (4.43) 


(4.44) 


(4.45) 






dll 


(4.46) 


4.9.3 Jet Mixing Model: 

There are many jet mixing models that can be used. In this program following models has 
been used. 

(a) Mass-Loss Model: 


In this model a specified fraction of the initial mass flow in the jet is assumed to mix with 
the main stream for each unit of distance along the jet. 

The mass flow remaining in the jet, mj , is given by: 


mj=m, 


1- 


Cs 


.A 


(4.47) 


J 


The entrainment constant, C, may be assigned any value between zero (no mixing) and 
one (instantaneous mixing). 

The other jet flow quantities are calculated directly from the residual mass flow, with the 
jet velocity, density, and temperature assumed uniform; 

A= ^ (4.48) 

P^^-,0 cos^j 

Mj = mjUj „ cos ^z?j (4.49) 


hj = mj 




(4.50) 


(b) Equivalent-Entrainment Model: 

This model starts with the actual mass flow in the jet, mj,act , as estimated fi'om Equation 
(4.36) for penetration jets or (4.44) for wall jets. In these equations the integration of 
mass flow is carried across the jet to a value of the jet width, t] , large enough that 
virtually all of the jet is included. The excess of this mass flow over the initial flow 
through the hole, (mj,act-mh), is the mass flow of gases that have been entrained into the 
jet from the main stream. It is then assumed that a fraction C of this entrained mass flow 
represents the jet mass flow that has been mixed into the main stream. The residual mass 
flow in the jet is then given by; 

(4-51) 

where C= 0 , no mixing 

C=1 , complete mixing when the computed flow of 

the actual jet reaches twice the initial jet flow. 

C= large number, implies instantaneous mixing 
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4.10 COMBUSTION ANALYSIS OF THE FLAME TUBE 

The equation of continuity, momentum, energy, and state are solved for a section of the 
flame tube between adjacent calculation Stations 1 and 2. Following sketch shows the 
flow contributions: 



Figure 4.7 Flame tube section 


4.10.1 Continuity Equation: 

The equation expressing the conservation of mass in the flame tube is: 

Pft.iUft.1 Aft,,- + 2 ™j.i+ 

residual jets residual jets entering jets 

Aft,2 “ ]^Aj,2 + ^^j,2 (4-52) 

all jets J all jets 

where min, = fuel burning rate 

4.10.2 Momentum Equation: 

The flame tube equations include terms for the entering and residual jets. Since the gas 
velocity in the flame tube is quite low, the friction loss has been ignored. 

The equation expressing the conservation of momentuin within a control volume 
bounded by the flame tube wall and by transverse planes at Station 1 and 2 is then: 

_ "f dA , 

P«,.Aft,,+mft,,Uft,,+ 2 , + 

entering jets residual jets 

= Pft.2Aft,2+mft.2Uft,2+ 

all jets 

where Mj = axial component of jet momentum 

= mjUj cos^2>j 

(p . = local angle between jet axis and flame tube center line 



64 


mft - mass flow of hot gases in flame tube. (This excludes 
unmixed air in the residual jets but includes the mass of fuel burned 
land 2) 


between Stations 


The integral pressure term is approximated by: 

dx ~~ 2 ~ ■^ft.iXpft.i + Pft, 2 ) 

The final form of flame tube momentum equation is then: 




+ E ^j.o+ E “^j,i 

entering jets residual jets 




•fl,2 


Kih. 





- Pft,2Aft.2 + mg_2Uft_2 + 2 (4.53) 

all jets 


4.10.3 Energy Equation: 

The equation expressing conservation of energy between adjacent calculation stations in 
the flame tube is written as follows, 




1 2 


m^,+ 


E hj.o«ij.o+ E 

entering jets residual jets 


j, + Jqdx 




m 


ft.2 + E^i 

all jets 


2"^j.2 


where the enthalpy of hot gasses is defined by: 


(4.54) 


h(T,)= |c,(T)dT 

Tb 

The specific heat Cp is correlated as a function of temperature by Equation (4.11); the 
enthalpy of jets, hj , is obtained from the jet mixing calculation. Equations (4.38), (4.46); 

and the heat release rate, q , from Equation (4.1 1). 

Heat loss from the hot gases to the flame tube walls has been ignored in Equation (4.49) 
to avoid the complications in the calculation procedure. 

The jet areas and mass flow rates are obtained from the results of jet mixing calculations. 
Equations (4.34), (4.36), (4.43), (4.44), (4.47), (4.48), and (4.51). 
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4.10.4 Equation of State: 

The equation of state may be written as follows: 

PiuTf ^ •' 

R= gas constant for air 

4.10.5 Method of Solution 

The mass flow at Station 2 is known from jet mixing results: 

mft,2 +E%c +Z“j.i “Z%2 

Since Equation (4.54) can not be solved explicitly for Tft ,2 , an iterative technique is used 
to arrive at a solution. In this equation the velocity terms make a small contribution to the 
total enthalpy. For the first stage in the iterative process, therefore, Ufi ,2 is set equal to Uft,i . 
In addition, the specific heat and heat release rate are computed for Tft ,2 =Tft,i. With these 
simplifications Equation (4.54) can be solved explicitly for a first estimate of Tfi^ . New 
values of the specific heat and heat release rate are then computed for the new Tft ,2 and 
Equation (4.54) is solved again. The process is repeated until Tft ,2 converges to an 
acceptable accuracy. 

Using the value of Tft ,2 so obtained equations are solved to obtain a quadratic equation in 
Pft,2 : 


Pft,2 




ft .2 




(4.57) 


where S^j .2 

the resulting value of pfi ,2 is substituted in the following continuity equation to obtain Uft ,2 : 

this value of Uft 2 is substituted in Equation (4.54) and a new value of Tft ,2 is found. This 
process is repeated until both Tfi ,2 and Uft ,2 have converged to within specified tolerances. 


4.1 1 LIMITATIONS OF TBOE AIR-FLOW SUBPROGRAM 


The more important limitations of the air-flow subprogram are as follows: 

1) The treatment of the primary zone as stirred reactor. 

2) The program does not calculate the fuel burning rate directly, but receives it as 
input and only serves to keep it below stoichiometric at any location. 

3) The jet mixing models suffer from a lack of experimental data upon which the 
choice of models and entrainment constants are based; in addition, the program 
does not relate the mixing rate of wall jets to the film cooling analysis carried out 
in the heat transfer section. 

4. 12 STRUCTURE OF THE SUBPROGRAM 

The air-flow subprogram consists of: 

1) A subroutine which steps up starting conditions and directs the air flow 
calculation through the other subroutines in the appropriate sequence. 

2) Subroutines for solving the equations for flow in the annulus and flame tube. 

3) Subroutines for calculating the heat addition due to fuel burning and the 
temperature in the primary zone. 

4) A subroutine for calculating the rate of mixing of the entering jets. 

The subroutines and flowcharts of air-flow subprogram are given in the Appendix A3. 



CHAPTETR 5 

COMBUSTOR HEAT- TRANSFER ANALYSIS 


5.1 INTRODUCTION 

The basic functions of the combustor dome and liners are to contain and guide the 
burning fiiel-air mixture from the point of fiiel injection to the turbine nozzle inlet, to 
meter airflow to the various zones within the combustor chamber and to protect engine 
structural components from the burning gases. Protection of the combustor dome and 
liners from the hot gases within the combustor is extremely important and challenging. 
The structure must be protected using cooling air at combustor inlet temperature. 

^2 LINER COOLING APPROACHES 

The established method for protecting combustor components from the hot products of 
combustion is to use film cooling, wWe a film of air at combustor inlet temperature is 
injected along the liner surface. Several different types of liner construction have been 
used for film-cooled liners as described below: 

5.2.1 Louvred liners 

This method of liner construction uses a simple sheet metal wall with a series of called 
louvres installed at intervals along the surface (Figure 5.1). 



iuPINGEMENT/FILM cooled 



transpiration cooled UNEH gipg 


hot side 

Figure 5.1 Liner cooling techniques 
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Liners of this type are lightweight, inexpensive and easily fabricated. 

Drawbacks associated with the louvred liner construction include short life due to 
cracking and poor control of cooling air flow. The base of the louvres is a prime location 
for crack initiation due to very high stress concentration at the point where the liner has 
been sheared during louvre installation. 

5.2.2 Cooling slots 

A majority of current combustors designs use some type of continuous circumferential 
annular slots spaced axially at intervals along the liner to inject a cooling film along the 
inner surface of the liner. Seven different slot configurations are shown in Figure 5.1 
below. 

5.3 ONE-DIMENSIONAL ANALYSIS 

Once the internal flow conditions are known, cooling flows are selected, and detailed 
geometry is specified, one-dimensional heat transfer analysis can be conducted at several 
different axial locations to determine whether peak liner temperatures are acceptable. If 
not, cooling flow and liner geometry are adjusted, and analysis is repeated until a 
satisfactory temperature profile is obtained. 

5.4 ANALYTICAL METHODS AND PROGRAM DEVELOPMENT 

The principal operations in the heat-transfer subprogram comprise the formation of the 
various heat flux components at the flame tube walls in terms of the wall temperature, the 
solution of the heat balance equation to determine the wall temperature and the 
evaluation of the wall heat fluxes. The heat transfer subprogram receives as input the 
geometry of the flame tube and casing fi’om the control subprogram and the axial 
distributions of velocity and temperature of the flame tube gases from the air flow 
subprogram. 

The primary objective of the heat transfer subprogram is to establish the axial 
distribution of temperature along the flame tube walls for various program options. 

5.4.1 Assumptions Made in Heat-Transfer Analysis 

^ Conditions are assumed to be constant in the circumferential direction for the 
annular geometry. 

> Steady-state conditions prevail. 

> The quantities supplied by the air flow subprogram are reasonable approximations 
to real conditions. 

> Flame and walls are gray radiators, i.e., spectral effects can be ignored. 

> Wall emissivities and absorptivities are constant. 

> Effect of wall to bulk temperature ratio can be neglected; fluid properties used in 
the heat transfer relations are evaluated at the bulk fluid temperature. 
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^ One dimensional radiation option assumes that the temperature gradients in the 
flame tube are small. 

> The temperature gradients along the walls are small and that the two walls are at 
approximately the same temperature. 

^ Unless the outer casing temperature distribution is specified in the program input, 
it is assumed that the casing temperature is close to the compressor discharge 
temperature. 

5.4.2 Heat Balance in Flame Tube 

The heat flux balance is performed on an elemental length of the flame tube wall. 
Conditions in the circumferential direction are assumed to be uniform. The heat flux 
components considered in the heat balance are: 

a) Radiation from the flame to the flame tube wall, Ri. 

b) Convection in the flame tube gases at the wall, Ci. 

c) Radiation from the flame tube wall to the outer casing, R 2 . 

d) Convection in the annulus air at the outside surface of the flame tube wall, C 2 

e) Longitudinal conduction in the flame tube wall, AK 

f) Radiation interchange between the flame tube walls, R3. 

Equating the heat gained by the element of wall to the heat lost, the following equation 
can be written; 


R, +Ci+AK = R3+C2+R3 (5.1) 

5.4.3 Radiation from the Flame to the Flame-Tube Wall (1-D Model) 

The one-dimensional radiation model describes the case where each element of wall 
receives radiation only from element of flame at the same axial position. This model is 
accurate when axial temperature gradients and end effects are small and the length-to- 
diameter ratio of the flame tube is large. 

The usual expression for radiant heat transfer between a gray flame at temperature 
Tft and a black container at temperature at Tw is; 

R,=4,T;-a.T:) (5.2) 

where, o = Stefan-Boltzmann constant 

Sft = emissisvity of flame at flame temperature 
ttft = absorptivity of flame to radiation at temperature Tw 
Sine the wall is not black, but has an absorptivity aw which will normally be 
between 0.8 and 1.0 a correction should be made for absorption, reflection, re-absorption 
and so on. McAdmas [31] suggests a factor (l-aw)/2 to allow for this effect. 

A further simplification is suggested by Lefebvre and Herbert [4]: 
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Then 


R.=54 + «kT.‘’(T,“-T“) 


5.4.4 Volume of Flame Element 


(5.3) 


The elements of flame considered have element length in the axial direction of AXft and 
inner and outer (flame-tube-wall) radii of T2 and ri respectively. Thus, the volume of each 
flame element is: 


AX, -4,^ 



(5.4) 


5.4.5 Emission per Unit Volume 


The radiation emitted by unit volume of flame temperature T* is, Hottel [31]: 




UPg^ 


-(5.5) 


Polb=0 


where, pQ = partial pressure of radiating gas (i.e., water vapor and carbon dioxide), 

ds. 


The quantity 




is mainly a function of temperature. For a typical value of 


b /Poli,=0 


pg/pft = 0.05 and for a non-luminous flame it can be approximated by [31], 

3.6 




(5.6) 


5.4.6 View Factor 


The sketch below illustrates the nomenclature of the axially symmetric configuration 
considered. 
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Figure 5.2 Flame tube view factor nomenclature 

The following equation was obtained for the view factor of the flame element radiating to 
an element of length AXw [28], 
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1 AX„ 

I— 


1 AX, 

1- ^ 

4(r3-rJ 

l/x“+fe 


4(r3-ri) 

Vx=+(r,-r,)\ 


(5.7) 


5.4.7 Transmittance 

The transmittance accounts for the radiative energy which is not absorbed by the 
intervening gas. It is a function of the distance the radiation travels through the gas, the 
gas composition and the gas temperature. The expression for transmittance, [18], is: 

„ 14.82 

PqS +14.82 

where s is the distance from element at considered point on the wall to element of the 
flame. 

Mean value of the transmittance for each axial location of the flame element to an 
outer wall element is given by: 


1 20 

■■—y 

20^ 


14.82 


14.82 + p<,JxH 


^3 h 


20 40 f' 


(5.9) 


5.4.8 Total Radiation 

Introduction of wall temperature, Tw, and effective wall absorptivity, (l + a^)/ 2, the 
radiation absorbed by the wall element from the flame element is given by: 


4ar„ 


(l + “w). 


^ de ^ 


dpoL 


AX,ii(r.'-r^)P™fe-T:) 


(5.10) 


7o 


Hence, the total radiation per unit time at a point on the outer flame-tube wall receiving 
radiation from all positions of the flame is: 

AX, (r.^-r,%(T;-T:) (5.11) 


2 all flamel 

pcsitions 


.dPcLJo 


The radiative heat flux per unit time and area of flame tube is therefore, given by: 


_ (1 + )a y, 

1 A.^ T% /mU 


AX ^ R. 3 all flame 
positions 


mean Jr G 


vdpoLy 


(5.12) 


A similar expression applies for radiation to the inner flame-tube wall. 


5.4.9 Radiation from Flame-Tube Wall to the Outer Casing 

The radiant heat flux from flame-tube wall to the outer casing is given by equation: 
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R>=F»°-(t:-T‘) (5.13) 

where Tc = temperature of casing. 

Fwc ~ overall interchange factor. 

For the case where only two surfaces are involved and where both surfaces have high 
emissivities, the overall interchange factor is approximately given by, [31]; 



(5.14) 


where, F^^ = black-surface overall interchange factor 

£c = emissivity of casing (assumed equal to absorptivity) 

Aw = surface area of flame-tube wall 
Ac = surface area of casing 

For radiation across an annular space, the black-surface exchange factor can be assumed 
to be equal to unity. With this simplification Equation (5.14) becomes. 

F„= (5.15) 

K 

Thus, from Equations (5.13) and (5.15) the radiation heat flux from the wall to the outer 
casing is: 

R,=:<t (t;-T‘) (5.16) 

- 

5.4.10 Radiation Exchange Between Flame-Tube Walls 

The net rate of heat transfer to a section of the flame-tube walls from all other ‘Visible” 
parts of the flame-tube wall is derived using the nomenclature shown in the sketch: 



Figure 5.3 View factor nomenclatures 
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Radiation emitted by an element 5xb which is received by element 6xa (neglecting any 
reflections) is: 

0£=2»i5x.F^.r(Ti-T:) (5.17) 

where, Fba = View factor for radiation from b to a 
5x3(rb-rJ^ 

^((rb-ra)' +(Xa-Xbyr 

r = transmittance between b and a 
14.82 


+k -hf +14. 
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(5.18) 


Radiation exchange per unit time and area of inner wall is given by: 

_ ost r.bXbfa-rJ^rfa-T^) 

I'a'tjrb +(xa -Xb)"!'" 

Thus the net rate of heat transfer per unit areas to the point on the inner wall under 
consideration, from all parts of the outer wall is: 


AR3 = 


(5.19) 


R, 


cs 


r S 


all b-wall 
positions 


r.tb-r.y+i*. 


(5.20) 


A similar expression applies for radiation from all parts of the inner wall to a point on the 
outer wall. 

5.4.11 Convection at the Inner Surface of the Flame-Tube Wall 

The correlation employed for the internal convective heat transfer at the flame-tube wall 
for subsonic flow through smooth tubes at high surface and fluid temperatures is [29], 


(5.21) 


C, = 0.023:^Re°*Pr°"(Tft - T„) 

where K = thermal conductivity of gas 

Dft = hydraulic diameter of the flame tube, equal to twice the distance between 
flame-tube walls. 


5.4.12 Entry-Length Effect 

When starting conditions on the wall are uniform, there is a developing boundary layer 
and the expression for the Nusselt number is usually multiplied by a factor for the entry- 
length effect. This is given as [31]: 


1 + 


Di 


-0.7 


(5.22) 
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5.4.13 Internal Convection with Film Cooling 

The contribution of internal convection to the wall heat balance is computed from 
Equation (5.21) only when film cooling and transpiration cooling are absent. For the case 
in which film cooling is employed. Equation (5.21) is modified by replacing the hot-gas 
temperature. Tfi by an adiabatic wall temperature, Thus: 


C,=0.023;^ReO*Pr°-^(T^-Tj 


D 


(5.23) 


The adiabatic wall temperature can be obtained from the following equations [8]: 

(5.24) 


T = 


T«-T 


X = 0.91 


U, 


cf 

s 0.8 




“+1.41 

X 



1 

1 


nO.5 


(5.25) 


where, Tcf = inlet temperature of cooling film 
Ycf = slot Avidth 

X = downstream distance of wall element from the cooling slot. 

5.4.14 Convection at the Outer Surface of the Flame-Tube Wall 

External convection from the flame-tube walls to the annulus air is governed by equation 
similar to Equation (5.21): 

C,=0.023iRe*‘Pr»MT.-T„) (5.26) 


5.4.15 Longitudinal Conduction Along flame-Tube Wall 


The longitudinal conduction effect can be estimated on the basis of average temperature 
gradients. For the inner flame-tube wall, the rate of heat flow to the section of wall (n) 
under consideration to an adjacent section (n+1) is given by: 
k^t„C^(T„„^, -T^ J 

H 


-X. 


(5.27) 


where kw = thermal conductivity of wall material 
tw = wall thickness 
Per unit area the heaf gain is: 


K = 


Xn+1-X„ 


x„^i-x„_, 


Combining this with a corresponding quantity for heat flow from the (n-1) element, the 
net heat received becomes: 


AK = 


"■n+l 


“X 


n-1 


T -T 

w.n+1 w.n 
Xn^l-Xn-l 


T -T 

W.n w.n-1 


x„ -X 


n-1 


(5.28) 



5.5 SOLUTION OF THE HEAT-BALANCE EQUATION 


When all of the heat-flux terms are considered, the heat-balance equation for an element 
of the flame-tube wall is: 

Ri+C, +AK=R2+C2+R3 (5,29) 

The longitudinal conduction term, AK, and the radiation exchange term, R3, are 
functions of the wall temperature distribution (rather than the local temperature of the 
wall element) and, therefore, an iterative method of solution of Equation (5.29) is 
required when these terms are considered. 

5.5.1 Non-Iterative Solution 

If we don’t consider AK and R 3 then a non-iterative solution of Equation (5.29) is 
possible. In this case Equation (5.29) becomes: 

R, +C, =R 2 +C 2 (5.30) 

Substituting for Ri from Equation (5.2) or (5.12); Cl from Equation (5.21), (5.22); for R 2 
from Equation (5.16); for Equation (5.25); then equation (5.30) becomes: 

DX+D 2 T.^'+D 3 T„=D, (5.31) 

where, Di, Da ,D3 and D4 are known. Newton’s approximation is used to solve Equation 
(5.31). 

For use in Newton’s approximation, a first estimate of Tw is obtained from one of 
the following equations: 

D = D 

T = T 

The chosen Tw is that which leads to the smallest value of s when Equation (5.31) is 
written in the form: 

DX+'DX' (5-32) 

With this initial estimate of Tw, Newton’s approximation is used to find a value of Tw that 
satisfies Equation (5.31). 

5.5.2 Iterative Solution 

The iterative heat-balance solution is only called for after a non-iterative calculation; an 
approximate wall-temperature distribution will therefore always be available. The 
iterative heat balance is repeated until the maximum temperature change between any 
two iteration cycles is less than specified tolerance. 



5.6 


LIMITATIONS OF THE HEAT TRANSFER SUBPROGRAM 


There are several significant limitations of the heat transfer subprogram: 

> The film cooling model is idealistic; there is a serious lack of experimental data 
upon which to base choice of the characteristic distance Xo used in the 
correlation; in addition, as mentioned above, the film cooling and wall jet mixing 
models are not related. 

> Data are lacking on the emissisvity of luminous flame. 

> No heat balance is performed to calculate the temperatures of the inner and outer 
casing walls. 

5.7 STRUCTURE OF THE SUBPROGRAM 

The heat-transfer subprogram consists of; 

1 . A subroutine for solving the non-iterative heat balance. 

2. A subroutine for solving the iterative heat balance. 

3. A group of library subroutines which provide gas properties, correlation for wall 
cooling and flame emissivity, and the solution of the fourth-power heat-balance 
equation. 

The subroutines and flowcharts of heat-transfer subprogram are given in the Appendix 

A4. 



CHAPTER 6 


RESULTS AND DISCUSSION 


6.1 INTRODUCTION 

This chapter presents the results obtained from the computer simulation of a gas turbine 
combustor. These results were the outcome of the techniques and methodology used in 
the previous chapters to analyze the combustor for flow and performance parameters. The 
results were compared with the available experimental data and 3D analysis for annular 
combustors of a modem fighter aircraft engine (Kaveri engine for LCA). 

There are many design aspects of a modem combustor and effort was to 
incorporate all these for the satisfactory operation of the engine, but this was not possible 
due to the complex nature of flow and heat transfer in the combustor. So a compromise 
was made between them. As for as 1-D analysis is concerned, the main parameters of 
interest are the pressure loss, the flow distribution and the heat transfer. As described in 
previous chapters, the analysis was carried out in three major sets, namely diffuser, main 
combustor, and heat transfer analysis. The results of the computer simulation program for 
various test mns and their comparison with the experimental data and 3D analysis 
illustrating aspects of the program are presented in the following three sections. 

6.2 DIFFUSER ANALYSIS RESULTS 

The diffuser analysis is mainly concerned with the checking of the validity of the 
analytical “stream tube” method. The results obtained in this subprogram were compared 
with the experimental data and 3D analysis on annular diffusers. 

Two types of combustor diffusers are taken into consideration, namely, K1 and 
K4 for the validation of the present diffuser subprogram. The analysis was carried out for 
two conditions, first for known mass flow split and second for unknown mass flow split. 
The inlet conditions and results are as follows: 

6.2.1 K1 -Combustor: 

Table 6.1 -Inlet Conditions 


Total Pressure (MPa) 

2.17 

MmsBssmmmm 


Air mass flow rate (kg/s) 

63.34 

Mach No. 

■HSESlHi 

Length (pre-diffuser) (m) 


Critical Shape Factor 

2.28 

Included Angle (degrees) 

6 

Mass flow inner annulus (%) 

42 

Mass flow outer annulus (%) 

45 

Mass flow core (%) 

13 







Table 6. 2 -Results of Pressure Loss and Mach number 


Pressure Loss 

1-D 

Known flow Split (13% core) 

0.14 

Unknovioi flow split 

0.15 

3-D 


0.18 

Mach No. 

1-D 

Known flow Split (13% core) 

0.2 

Unknown flow split 

0.2 

3-D 


0.17 

Ideal Pressure ] 

Recovery Coefficient (Cpi ) 

0.971 

Pressure Recovery Coefficient (Cp ) 

0.645 

Diffuser Effectiveness (^) 

0.656 


6,2.2 K4 -Combustor; 


Table 6.3 - Inlet Conditions 


Total Pressure (MPa) 

2.17 

Total Temperature (K) 

756.50 

Air mass flow rate (kg/s) 

63.34 

Mach No. 

0.34 

Length (pre-diffuser) (m) 

0.065 

Critical Shape Factor 

2.28 

Included Angle (degrees) 

7 

Mass flow inner annulus (%) 

42 

Mass flow out^r annulus (%) 

45 

Mass flow core (%) 

13 
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Table 6.4 - Results of Pressure Loss and Mach number 


Pressure Loss 

1-D 

Known flow Split (13% core) 

0.17 

Unknown flow split 

0.18 

3-D 


0.2 

Mach No. 

1-D 

Known flow Split (13% core) 

0.23 

Unknown flow split 

0.23 

3-D 

1 


0.216 

Ideal Pressure 

Recovery Coefficient (Cpi ) 

0.972 

Pressure Recovery Coefficient ( Cp ) i 

0.945 

Diffuser Effectiveness (^) 

0.983 


Tables 6.2 and 6.4 show the pressure loss, outlet Mach number, pressure recovery 
coeflficient, and effectiveness of the diffuser as obtained from diffuser calculations. These 
results are in good agreement with the 3D analysis, which had been verified and validated 
with the experimental data for the same diffuser under same inlet conditions. The 
differences among, 3D analysis, experimental data and ID is due to the fact that when the 
stream tube method is used between the annuli, the calculated boundary layer blockage 
decreases considerably between the outlet of pre-diffiiser and inlet to the aimulus which 
results in lower mixing losses. 

The pressure losses in the diffusing passages influence the mass flow split at the 
snout or cowl. The flow split at the snout in turn influences the temperature level in the 
flame tube. The diffuser analysis has been carried out by the ‘stream tube method’ 
assuming a uniform inlet velocity profile; to account for the fact that the actual velocity at 
the wall was about 20 per cent greater than that in the mid stream, the test cases were run 
with a blockage of 0.8, and a corresponding shape factor of 1.24. Since the deceleration 
of the flow between the diffuser walls is greater for the non-uniform velocity profile, a 
larger boundary layer blockage should be used with a uniform velocity profile to simulate 
this effect. The only effect which the computations apparently do not accurately predict is 
the effect of high Mach numbers on the pressure-loss factor. 

6.3 COMBUSTOR FLOW ANALYSIS RESULTS 

The combustor flow analysis is mainly concerned with the analysis for air mass flow split, 
inner and outer annulus pressure losses, flow distribution through various cooling holes 
and slots and outlet Mach numbers in annuli and flame tube. The inlet conditions are the 
same as given in the Table 6.1 for K1 combustor and Table 6.3 for K4 combustor. The 
results pertaining to the combustor flow for K1 and K4 type of combustors are given 

below. 
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6.3.1 Kl- Combustor: (same inlet conditioned as mentioned in Table 6.1) 


Table 6.5 — Results of Liner Mass Flow Distribution 



1-D 

3-D 

Known 

Split(13% 

Core) 

Unknown 

Split 

Outer Annulus 

% 

% 

% 

Primary Zone 

10.92 

10.97 

10.9 

Cooling Ring No. 3 

2.71 

2.72 

2.9 

Secondary Zone 

2.43 

2.44 

3.2 

Cooling Ring No. 5 

2.68 

2.68 

2.7 

Dilution Zone 

8.31 

8.34 

10.5 

Cooling Ring No. 7 

2.91 

2.92 

2.9 

Cooling Ring No.9 

2.85 

2.86 

2.6 

Cooling Ring No.l 1 

2.84 

2.85 

2.9 

Inner Annulus 




Primary Zone 

11.34 

11.28 

12.5 

Cooling Ring No.4 

1.92 

1.92 

1.6 

Secondary Zone 

2.59 

2.59 

3.0 

Cooling Ring No.6 

1.94 

1.94 ^ 

1.5 

Dilution Zone 

10.04 

10.02 

10.1 

Cooling Ring No. 8 

2.02 

2.02 

2.2 

Cooling Ring No. 10 

2.28 

2.28 

2.0 

Cooling Ring No. 12 

2.45 

2.46 

2.2 
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Table 6.6 - Results of Pressure Loss and Mach number 



Inner 

Annulus 

Outer 

Annulus 

Overall 
(Flame 
Tube Exit) 

Pressure 

Loss 

1-D 

Known Split 
(13%Core) 

2.04 

2.7 

6.8 

Unknown 

2.05 

2.7 

6.86 

3-D 


1.29 

1.97 

6.83 

Mach No 

1-D 

Known Split 
(13%Core) 

0.111 

0.096 

0.124 

Unknown 

0.112 

0.095 

0.123 

3-D 


0.105 

0.0984 

0.119 


Table 6.7 - Results of Mass Flow Split 



Inner Annulus 
(%) 

Outer Annulus 
(%) 

Core 

(%) 

1-D 

Known Split 
(13%Core) 

41 

44 

15 

Unknown 

41.4 

43.6 

15 

3-D 

42.3 

44.7 

13 





6.3.2 K4- Combustor : (same inlet conditions as given in Table 6.3) 

Table 6.8 - Results of Liner Mass Flow Distribution 



1-D 

Experimental 

3-D 

Known 

Split(13% 

Core) 

Unknown 

Split 

Outer Annulus 

% 

% 

% 

% 

Primary Zone 

11.00 

10.840 

12.40 

10.7 

Cooling Ring No. 3 

2.790 

2.760 

2.764 

3.1 

Secondary Zone 

2.434 

2.400 

2.450 

3.2 

Cooling Ring No. 5 

2.764 

2.740 

2.764 

2.7 

Dilution Zone 

8.474 

8.380 

8.840 

10.6 

Cooling Ring No. 7 

3.016 

3.000 

2.764 

2.9 

Cooling Ring No.9 

2.905 

2.880 

2:764 

2.7 

Cooling Ring No. 1 1 

2.856 

2.830 

2.764 

2.9 

Inner Annulus 





Primary Zone 

12.030 

11.910 

14.640 

12.4 

Cooling Ring No.4 

1.923 

1.910 

2.234 

1.6 

Secondary Zone 

2.342 

2.320 

2.900 

3.0 

Cooling Ring No. 6 

1.923 

1.910 

2.234 

1.5 

Dilution Zone 

9.463 

9.360 

10.18 

10.1 

Cooling Ring No. 8 

2.014 

2.000 

2.234 

2.2 

Cooling Ring No. 10 

2.167 

2.150 

2.234 

1.9 

Cooling Ring No . 1 2 

2.344 

. 

2.330 

2.234 

2.2 
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Table 6.9 - Results- Total Pressure Loss and Mach number 



Inner 

Annulus 

Outer 

Annulus 

Overall 
(Flame 
Tube Exit) 

Pressure 

Loss 

1-D 

Known Split 
(13%Core) 

2.14 

3.12 

8.44 

Unknown 

2.14 

3.12 

8.43 

3-D 


2.32 

2.71 

9.17 

Mach No 

1-D 

Known Split 
(13%Core) 

0.125 

0.108 

0.16 

Unknown 

0.125 

0.109 

0.16 

3-D 


0.1343 

0.118 

0.1373 


The results shown in Tables 5, 6 and 7 pertain to the K1 combustor. The results are in 
agreement with the 3D analysis and hence with experimental results since the 3D analysis 
have already been validated with the experimental results. There is a little difference in 
the dome mass flow split. .This variation may be due to the fact that ID analysis 
calculates the mass flow split based on the flow area. Since the discharge coefficient of 
holes have influence over the flow mass rate, its effect on flow split can not be 
overlooked. An accurate value of discharge coefficient for different types of holes under 
different flow conditions would surely increase the accuracy of the results obtained from 
the computer program. 

Table 6.8 and 6.9 show the results of K4 combustor. Table 6.8 compares ID and 
3D predictions with the well established experimental data. One can clearly see that 
thelD predictions are as close to the experimental results as the 3D predictions. 


6.4 HEATS-TRANSFER RESULTS 

There are many heat transfer parameters that affect the functioning of an aircraft engine 
combustor like flame tube wall temperature, film cooling effectiveness and the heat 
transfer to the annulus air mainly due to convection and radiation. The flame tube wall 
temperature is the most critical and important parameter for a combustor design. It 
imposes a certain limit on the maximum temperature rise of liner wall material to avoid 
distortion and to keep the mechanical strength adequate enough. 

1-D analysis has been carried out for various heat transfer and cooling parameters 
of the combustion chamber. Here only flame tube wall temperature results are presented 
due to its importance over other parameters. 



The predictions of the heat-transfer subprogram for K1 and K4 have been 
compared with data presented by NASA [35], The casing temperature was calculated 
assuming natural convection and radiation using room-temperature as the reference. (The 
actual heat transfer coefficients were obtained from McAdams [31]). Heat transfer from 
the casmg to the annulus air was ignored. The flame tube wall absorptivities and 
emissivities were taken as 0.85; the casing emissivity was taken as 0.8 

Figure 6.1 shows the temperature distribution of flame tube wall for a NASA 
combustor. The geometric configuration of this combustor is quite simple as compared 
to K1 and K4 combustors. Unlike K1 and K4, this combustor doesn’t have any dump 
region in the diffuser part. Its inlet conditions are also different. Still this result has been 
used for the validation of the computer program for heat transfer in the absence of other 
experimental data. This comparison, at least, gives a trend and range of predicted values 
from the present simulation. 



Axial Distance (mm) 

Figure 6.1 NASA-Liner wall temperature distributions. 



Figures 6.2 and 6.3 pertain to the K1 and K4 combustors, respectively, and give the linear 
temperature distribution obtained from ID prediction. In spite of different geometry and 
inlet conditions of NASA combustor, agreement of these results with Fig. 6.1(NASA) is 
reasonably good. The virtually instantaneous mixing of the penetration jets results in a 
sharply falling flame temperature immediately downstream of each hole row. As the wall 
jets don’t mix at all, the amount of gas available in the tube is reduced by the contents of 
the wall jets and the resulting spatial variation of flame-temperature is considerably 
higher than it would have been, had the wall jets mixed with the main stream in the flame 
tube. 



Figure 6.2 Kl-Liner wall temperature distributions. 

As evident from the plots, the temperatures of the flame tube walls follow the 
general variation of the flame temperature reasonably closely up to an axial position of 
about 370 mm. as it would be expected with a 1-D calculation of radiation from the flame 
tube walls. Downstream of this position, the sharply dropping radiative heat transfer from 
the flame is roughly balanced by a drop in convective heat transfer in the annuli caused 


by the sharply reduced air velocities in the annuli downstream of the dilution holes. Near 
the end of the combustor the air velocities in the annuli fall nearly to zero, causing the 
wall temperature to rise sharply. 

A closer look on the curves of Fig. 6.2 and 6.3 shows certain variations of the 
predicted values with NASA results of Fig. 6.1. The inconsistency may be due to; 

1 The difference in geometry and inlet condition although the basic shape is same. 

2. Simpler assumptions about the dimensions of the test combustor. 

3 . Simpler assumptions about the heat loss from the casing. 

4. Two-dimensional effects which become increasingly important as the convective 
heat transfer components are reduced. 

The casing temperature exerts an important influence on the flame tube wall temperature, 
since radiation froih the flame tube wall to the casing becomes significant in the absence 
of film cooling. 



Axial Distance (mm) 

Figure 6.3 K4-Liner wall temperature distributions. 




The most striking differences between the above curves are in their slopes. There 
appears to be a dominant two-dimensional heat transfer mode that is causing the 
experimentally measured wall temperatures to fall off rapidly with distance downstream. 
This effect is most apparent in the radiation-only results. As the amount of fuel available 
for burning in the primary zone does not vary from case to case, and the total temperature 
of the entering air is almost constant, the effect of a decrease in the flow through the 
dome into the primary zone is to raise the temperature of the combustion products. From 
the results of the present computer calculations, it is concluded that the film and 
convective cooling models used are producing results that agree with the NASA values. 


6.5 CONCLUSIONS 

Operation of the complete computer program has been checked by running two test cases 
having different geometry, but identical inlet conditions. The results of the diffuser, 
airflow and heat transfer test cases have been discussed in their respective sections 
separately. The comparison of results with well established data reveals the accuracy and 
validity of the computer program for ID analysis of a modern annular combustor. 

The attempt was to predict the maximum number of flow and performance 
parameters for an annular combustor using ID analysis for reasonably accurate results. 
The results obtained from analysis are presented graphically and in Tabular forms. The 
main conclusions of this work are as follows; 

1. A computer program for analyzing the airflow and heat transfer characteristics of 
ann ular gas turbine combustors has been developed. 

2. The diffuser, air-flow, and heat-transfer sections of the computer program were tested 
individually by comparing their predictions for appropriate quantities with available 
experimental data and 3D analysis in the open literature. The correlation between 
calculated and experimental values is judged to be reasonably good. 

3. The complete computer program was tested for satisfactory operation by running it 
for two test cases, using different geometry, but same inlet conditions and program 
options as specified by GTRE. The main trends of the results, as inferred by 
comparing changes in pressures, flows, and temperature distributions from case to 
case, are also in agreement with the 3D analysis as well as with the available 
experimental data. 



6.8 SUGGESTIONS FOR FUTURE WORK 


Present computer modeling can be made more powerful and accurate by introducing the 
analysis for those regions in which sudden expansion or eddy losses are predominant. 
Dump diffuser region analysis comes in the same category. As ID model is not capable 
of capturing the eddy and turbulence losses, which are at least 2D in nature; therefore ID 
computer simulation could be coupled with 2D or 3D modeling specifically for the dump 
region. This would help in providing more realistic input conditions for the air-flow and 
heat-transfer subprograms. 

The results of the stream tube diffuser calculation are quite good in view of the 
gross assumptions made. If the equations are examined further and the potentialities and 
limitations of the method explored, this technique could prove to be a powerful new tool 
for more general applications. 

It would be useful to make a more exhaustive survey of the influence of certain 
program input parameters. For example, critical shape factor, boundary layer blockage, 
entrainment function and constant, hole discharge coefficient, etc, are few parameters to 
which program results are very sensitive. By experiment, a “most reasonable” choice of 
these parameters may be found and accuracy of the simulation results can be increased 
significantly. 
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APPENDIX 


SUBROUTINES AND FLOWCHARTS 


A1 INTRODUCTION 


This section is meant to show how the modeling methods described in Chapter 2 through 
Chapter 5 are fitted into the computer program. As already explained that the computer 
modeling program consist of three independent subprograms, viz., diffuser, airflow and heat 
transfer subprograms, each of which contains a family of the subprograms 



Figure A1 Flowchart for Main Program 


Figure A1 shows the main program flowcharts. It calls subroutine DIFLOW for the diffuser 
portion of calculation under diffuser subprogram, subroutine AIRFLOW for the air-flow 
calculation under airflow subprogram and subroutine HEAT for heat transfer calculation of 
combustor under heat transfer subprogram. Figure A2 presents flowchart for overall program 
including as a whole including reading of input data and writing out it to make sure that the 
input data are read in correctly. Geometrical calculations of the combustor are also carried out 
to supply required geometric input data to various subroutines in appropriate form for further 
calculations in the subroutine. The subroutines and respective flowcharts are given in brief in 
the succeeding sections. 
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Figure A2 Overall Flowchart for Entire Program 


Heat-transfer Subnroaram 
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A2 DIFFUSER SUBPROGRAM 

The diflEliser subprogram calculates the flow properties in the diffuser from the 
compressor outlet to the first calculation point in each annulus. The subroutines and 
Flowcharts are in the following subsections. 

A2.1 Diffuser Subprogram Subroutines 

Various subroutines pertaining to this subprogram are as follows: 

Subroutine DIFLOW: This subroutine has three main functions- 

1 . To organize the calculation procedure used in the diffuser subprogram. 

2. To calculate the flow properties in the two annuli between the cowl and outer 
casing if the mixing equation is being used. 

3 . To calculate the flow properties beyond the outlet of pre-diffuser (a fudge 
factor has been used for the results obtained in this region as to compensate for 
the eddy and turbulence losses to make available more realistic input values for 
the next calculation stations). 

This subroutine has three entry points splitting the program into three separate parts. The three 
parts performs the following fiinction: 

(i) The first part organizes the calculation position between the pre-diffiiser inlet 
to the cowl inlet. A correction for the passage curvature is also calculated. 

(ii) The second part, following DIFLW, the calculation procedure between the 
dome and annulus diffusing passage is organized and flow properties are 
calculated at the outlet. 

(iii) The third part, following Entry DIFLW2, calculates the flow in the diffuser 
with a given mass flow split using the stream tube method throughout. 

Subroutine TUBCTS: This subroutine is called by DIFLOW. The subroutine does the 
following: 

1 . Calculates the inlet static pressure and the inlet stream tube properties. 

2. Calculates the reference velocity, Mach number, and dynamic head. 

3 . Sets up an iteration loop to calculate the flow properties and the boundary layer 
characteristics. 

4. Calculates the outlet velocity profile and the diffuser performance. 

Subroutine TUBSTA: It is called by DIFLOW. This subroutine is called before the iteration 
on mass flow commences, and its function is to set up geometric data in a form 
that can be used by subroutine TUBSAl. 

Subroutine TUBSAl: This subroutine is called by DIFLOW to perform a stream tube 
calculation in the diffuser part beyond pre-diffuser outlet to: 
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1 . Calculate the stream tube properties. 

2. Calculate the boundary layer characteristics. 

3 . Calculates the diffuser performance. 

Subroutine TUBFWl: It is called by subroutine DIFLOW during the iteration on the mass 
flow split. This calculates; 

1 . The geometric positions of the flow split 

2. The area mean velocity and weighted mean velocity for the flow streams into 
each annuls and the cowl or snout. 

Subroutine TUNEIN: This is called by TUBCTS and TUBSAl to calculate the input 

quantities for the stream tube calculation. In this the following calculations are 
performed- 

1 . If the static pressure is not given (i.e., if TUBEIN is called by TUBCTS), 
calculate the static pressure. 

2. Split the flow into a number of stream tubes and calculates the total pressure 
for each stream tube. 

3 . Calculates the flow area for each stream tube when the flow is accelerated 
isentropically to Mach. Number 1. 

Subroutine TUBANL: This is called by TUBCTS and TUBAS 1 to calculate- 

1 . The pressure at each calculation point in the diffuser, Mach number. 

2. The velocity and density at each calculation point. 

3. Outlet velocity profile 

Subroutine NEWRAD: This is called by TUBCTS and TUBSAl to perform- 

1 . Position of flow separation (depending upon critical shape factor supplied 
from outside to the program). 

2. Checks the solution for convergence. 

3 . Calculates a new guess for the boundary layer displacement thickness. 

A2.2 Flowchart for Diffuser Subprogram 

The overall flow chart for the diffuser subprogram is given in the Figure A3. 
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Figure A3 Overall Flowchart for Diffuser Subprogram 

















A3 AIR FLOW SUBPROGRAM 

The air-flow subprogram calculates the flow properties in the flame tube and annuli 

from the first calculation point, which corresponds the first hole row on the wall as distinct 

from the dome, to the end of the combustor. It also organizes the mass flow split iteration. 

A3.1 Air-Flow Subprogram Subroutines: Subroutines of this subprogram are 

Subroutine AIRFLO: This subroutine controls the air-flow calculations. It functions are- 

1 . An initial mass flow split is calculated based on the total hole area in the dome 
on each annulus. 

2. Calculates the initial conditions on the dome at the start of each annulus. 

3 . Calculates the primary zone pressure and temperature also initial conditions for 
the flame tube calculations. 

Subroutine DISJET: This subroutine determines jet discharge parameters which are; 

1. Initial jet-air mass flow rate, discharge jet angle, and discharge coefficient. 
This jet information is used by subroutine JETMIX. 

2. Initial jet-air mass flow rate, jet momentum, and axial jet momentum and jet 
enthalpy. This jet information is used by subroutine EQUAN and EQUFT. 

Subroutine PRTEMP; This subroutine does the following functions: 

1 . Calculates conditions within the primary zone by solution of a simplified form 
of the flame tube equations, assuming the primary zone acts as a stirred reactor. 
Flame tube conditions are thus set up for the secondary-hole calculation point, 
so that from then on the flame tube equations can be solved in the normal way. 

2. This also calculates the fraction of air entering the primary holes which 
recirculates upstream. 

Subroutine EQUAN: This subroutine solves the equations for flow in the annuli. The main 
equation is a quadratic in annulus velocity at the point of calculation. The rest 
of subroutine is concerned with evaluating the constants in this equation. 

Subroutine EQUFT: This subroutine solves the equations for flow in the flame tube down 
stream of the primary zone. The solution follows an iterative procedure. 
Calculations are performed for a control volume bounded by two adjacent 
calculation stations (say 1 and 2). 

Subroutine HEATAD: This subroutine provides EQUFT with the heat release terms in the 
energy equation. 

Subroutine JETMIX: This subroutine provides with the characteristics of the residual jets. 

Wall jets and penetration jets are calculated separately using the same model 
but different values of the mixing constants. 

A3.2 Flowchart for AirOflow Subprogram: The overall flowchart for 

diffuser subprogram is given in Figure A4. 
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Figure A4 Overall Flowchart for Air-Flow Subprogram 



















A4 HEAT TRANSFER SUBPROGRAM 

The heat transfer subprogram calculates the temperature distribution in the liner wall, 
film effectiveness and the rate of heat transfer among various components of the m ain 
combustor. Main subroutines and flowchart of the program is given in the following 
sections. 

A4.1 Heat Transfer Subprogram Subroutines: The main subroutines of the 
subprogram are as follows. 

Subroutine HEATl: This subroutine carries out the following- 

1 . This subroutine performs non-iterative heat transfer calculation. 

2. The coefficients in the heat-balance equation are evaluated, using 
subroutines EEFT, COOL, and PROP. 

Subroutine HEAT2: This subroutine does the following fonctions- 

1 . This subroutine carries out an iterative heat-transfer calculation. 

2. The subroutine also writes out the results for the heat transfer subprogram 
in the readable form. 

Subroutine TWALL: This subroutine uses Newton’s approximation to solve the heat 
balance equation, which is in the following form; 

DiTw'^+D2Tw"-^+D3Tw = D4 

Subroutine EEFT: This subroutine calculates the emissivity of the flame. 

Subroutine PROP: This subroutine obtains average values for the specific heat, thermal 
conductivity, and the dynamic viscosity of the gas mixture. 

Subroutine COOL; This subroutine calculates wall cooling parameters according to the 
film cooling technique. 


A4.2 Flowchart for Heat-Transfer Subprogram: 

The overall flow chart for the heat transfer subprogram is given in the Figure A5. 



yy 



Figure A5 Overall Flowchart for Heat-Transfer Subprogram 













